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Abstract 

Within thunderstorms electrons can gain energies of up to hundred(s) of 
MeV. These electrons can create X-rays and gamma-rays as Bremsstrahlung 
when they collide with air molecules. Here we calculate the distribution 
r-| \ of angles between incident electrons and emitted photons as a function of 

^ electron and photon energy. We derive these doubly differential cross-sections 

by integrating analytically over the triply differential cross-sections derived 
by Bethe and Heitler; this is appropriate for light atoms like nitrogen and 
oxygen and for electron energies between 1 keV and 1 GeV. We also discuss 
some simplifying limit cases, and we derive some simple approximation for 
the most probable scattering angle. 

We also provide cross sections for the production of electron positron 
pairs from energetic photons when they interact with air molecules. This 
process is related to the Bremsstrahlung process by some physical symmetry. 
Therefore the results above can be transferred to predictions on the angles 
between incident photon and emitted positron, again as a function of photon 
and positron energy. We present the distribution of angles and again a simple 
approximation for the most probable scattering angle. 

Our results are given as analytical expressions as well as in the form of a 
C++ code that can be directly be implemented into Monte Carlo codes. 

Keywords: Bremsstrahlung, pair production, analytical integration of 
Bethe-Heitler equation, distribution of scattering angles 
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1. Introduction 



1.1. Flashes of gamma-rays, electrons and positrons above thunderclouds 

Terrestrial gamma ray fiaslies (TGFs) were first observed above tliunder- 
clouds by tlie Burst and Transient Source Experiment (BATSE) (Fisliman 
et al., 1994). It was soon understood that these energetic photons were gen- 
erated by the Bremsstrahlung process when energetic electrons collide with 
air molecules ( iFishman et al.. 19941 iTorii et al.. 20041 ): these electrons were 
accelerated by some mechanism within the thunderstorm. Since then, mea- 
surements of TGF's were extended and largely refined by the Reuven Ramaty 
Energy Solar Spectroscopic Imager (RHESSI) (Cummer et al., 2005; Smith 
et al., 2005, Grefenstette et al., 2009, S mith et al., 2010: Lu et al., 2011), by 
the Fermi Gamma-ray Space Telescope ( Briggs et al.. 20l"ol ). by the Astror- 
ivelatore Gamma a Immagini Leggero (AGILE) satelli te which recently mea- 
sured TGFs with quantum energies of up to 100 MeV flMarisaldi et al. 2010l : 



Smith et al.. 20 111 : iTavani et al.. 201 ll ) , and by the Gamma- Ray Observati on 



of Winter Thunderclouds (GROWTH) experiment (ITsuchiya et al.. 20111 ). 

Hard radiation was also measured from approaching lightning leaders 
(Moore et al., 2001; Dwyer et al., 2005); and there are also a number of 
laboratory experiments where very energetic photons were generated during 
the streamer-leader stage of discharges in open air (Stankevich and Kalinin, 
1967; Dwyer et al., 2005b; Kostyrya et al., 2006; Dwyer et al., 2008b; Nguyen 
et al., 2008; Rahmen et al., 2008; Rep'ev and Repin, 2008; Nguyen et al., 
2010; March and Montanya, 2010; Shao et al., 2011). 

Next to gamma-ray flashe s, flashes of energeti c electrons have been de- 
tected above thunderstorms ( Dwyer et al.. 2008bl ): they are distinguished 
from gamma-ray flashes by their dispersion and their location relative to 
the cloud - as charged particles in sufficiently thin air follow the geomagnetic 
field lines. In December 2009 NASA's Fermi satelli te detected a subs tan- 
tial amount of positrons within these electron beams ( iBriggs et al.. 20111 ). It 
is now generally assumed that these positrons come from electron positron 
pairs that are generated when gamma-rays collide with air molecules. 

Two different mechanisms for creating large amounts of energetic elec- 
trons in thunderclouds are presently discussed in the literature. The older 
suggestion is a relativistic run-away process in a rather homogeneous electric 
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field inside the cloud (Wilson, 1925; Gurevich, 1961; Gurevich ct al., 1992; 
Gurcvicli, 2001 Dwyer, 2003, 2007; Milikh and Roussel-Dupre, 2010). 

More recently research focusses on electron acceleration in the streamer- 
leader process with its strong local field enhancement (Moss et al., 2006; Li 
et al., 2007; Chanrion and Neubert, 2008; Li et al, 2009; Carlson et al, 2010; 
Celestin and Pasko, 2011; Li et al, 2012). 

1.2. The need for doubly differential cross-sections 

Whatever the mechanism of electron acceleration in thunderstorms is, ul- 
timately one needs to calculate the energy spectrum and angular distribution 
of the emitted Bremsstrahlung photons. As the electrons at the source form a 
rather directed beam pointing against the direction of the local field, the elec- 
tron energy distribution together with the angles and energies of the emitted 
photons determine the photon energy spectrum measured by some remote 
detector. The energy resolved photon scattering angles are determined by 
so-called doubly differential cross-sections that resolve simultaneously energy 
hoj and scattering angle ©j of the photons for given energy Ei of the incident 
electrons. The data is required for scattering on the light elements nitrogen 
and oxygen with atomic numbers Z = 7 and Z = 8, while much research in 
the past has focussed on metals with large atomic numbers Z. The energy 
range up to 1 GeV is relevant for TGF's; we here will provide data vahd for 
energies above 1 keV. 

As illustrated by Fig. 1, the full scattering problem is characterized by 
three angles. The two additional angles 0/ and $ determine the direction 
of the scattered electron relative to the incident electron and the emitted 
photon. The full angular and energy dependence of this process is determined 
by so-called triply differential cross-sections. A main result of the present 
paper is the analytical integration over the angles G/ and $ to determine 
the doubly differential cross-sections relevant for TGF's. 

As the cross-sections for the production of electron positron pairs from 
photons in the field of some nucleus are related by some physical symmetry 
to the Bremsstrahlung process, we study these processes as well; we provide 
doubly differential cross-sections for scattering angle B_|_ and energy of 
the emitted positrons for given incident photon energy hw and atomic number 
Z. 
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Figure 1: Parametrization of the Bremsstrahlung process: Momenta of incident electron 
Pi, scattered electron p/ and emitted photon ?ik form the angles 0^ = <(pi,k) and 
0/ = <(p/, k), and $ is the angle between the planes spanned by the vector pairs (pi, k) 
and (p/,k). The scattering nucleus has atomic number Z. 

1.3. Available cross-sections for Bremsstrahlung 

Our present understanding of Bremsstrahlung and pair production was 
largely developed in the first half of the 20th century. It was first calculated 
by Bethe and Heitler (1934). Important older reviews are by Heitler (1944), 
by Hough (1948), and by Koch and Motz (1959). We also used some re- 
cent text books (Greiner and Reinhardt, 1995; Peskin and Schroeder, 1995); 
together with Heitler (1944) and Hough (1948), they provide a good intro- 
duction into the quantum field theoretical description of Bremsstrahlung and 
pair production. The calculation of these two processes is related through 
some physical symmetry as will be explained in chapter [3J 

As drawn in Fig. 1, when an electron scatters at a nucleus, a photon with 
frequency u can be emitted. The geometry of this process is described by the 
three angles Gj, 9/ and $. Cross sections can be total or differential. Total 
cross sections determine whether a collision takes place for given incident 
electron energy, singly differential cross sections give additional information 
on the photon energy or on the angle between incident electron and emitted 
photon, and doubly differential cross sections contain both. Triply differential 
cross sections additionally contain the angle at which the electron is scattered. 
As two angles are required to characterize the direction of the scattered 
electron, one could argue that this cross section should actually be called 
quadruply differential, but the standard terminology for the process is triply 
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differential. 

Koch and Motz (1959) review many different expressions for different 
limiting cases, but without derivations. Moreover, some experimental results 
are discussed and compared with the presented theory. Bethe and Heitler 
(1934), Heitler (1944), Hough (1948), Koch and Motz (1959), Peskin and 
Schroeder (1995), Greiner and Reinhardt (1955) use the Born approximation 
to derive and describe Bremsstrahlung and pair production cross sections. 

Several years later new ansatzes were made to describe Bremsstrahlung. 
Elwert and Haug (1969) use approximate Sommerfeld-Maue eigenf unctions 
to derive cross sections for Bremsstrahlung under the assumption of a pure 
Coulomb field. They derive a triply differential cross section and beyond 
that also numerically a doubly differential cross section. Furthermore they 
compare with results obtained by using the Born approximation. They show 
that there is a small discrepancy for high atomic numbers between the Bethe- 
Heitler theory and experimental data, and they provide a correcting factor 
to fit the Bethe-Heitler approximation better to experimental data for large 
Z . However, they only investigate properties of Bremsstrahlung for Z = 13 
(aluminum) and Z = 79 (gold). 

Tseng and Pratt (1971) and Fink and Pratt (1973) use exact numerical 
calculations using Coulomb screened potentials and Furry- Sommerfeld-Maue 
wave functions, respectively. They investigate Bremsstrahlung and pair pro- 
duction for Z = 13 and for Z = 79 and show that their results with more 
accurate wave functions do not fit with the Bethe Heitler cross section ex- 
actly. This is not surprising as the Bethe-Heitler approximation is developed 
for low atomic numbers Z and for Z dependent electron energies as discussed 
in section 12.21 

Shaffer et al. (1996) review the Bethe Heitler and the Elwert Haug the- 
ory. They discuss that the Bethe Heitler approach is good for small atomic 
numbers and give a limit of Z > 29 for experiments to deviate from theory. 
They calculate triply differential cross sections using partial-wave and mul- 
tipole expansions in a screened potential numerically for Z = 47 (silver) and 
Z = 79 and compare their results with experimental data. Actually their 
results are close to the Elwert Haug theory which fits the experimental data 
better than their theory. 

Shaffer and Pratt (1997) also discuss the theory of Elwert and Haug 
(1969) and compare it with the Bethe Heitler theory and, additionally, with 
the Bethe Heitler results multiplied with the Elwert factor and with the exact 
partial wave method. They show that all theories agree within a factor 10 
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in the keV energy range, and that the Elwert-Haug theory fits the exact 
partial wave method best. However, they only investigate Bremsstrahlung 
for atomic nuclei with Z = 47, 53 (iodine), 60 (neodymium), 68 (erbium) 
and 79, but not for small atomic numbers Z = 7 and 8 as relevant in air. In 
summary, Elwert and Haug (1969), Tseng and Pratt (1971), Fink and Pratt 
(1973), Shaffer et al. (1996) and Shaffer and Pratt (1997) calculate cross 
sections for Bremsstrahlung and pair production for atomic numbers Z = 13 
and Z > 47 numerically, but not analytically, and they do not provide any 
formula or data which can be used to simulate discharges in air. 

The EEDL database consists mainly of experimental data which have 
been adjusted to nuclear model calculations. For the low energy range Geant4 
takes over this data and gives a fit formula. The singly differential cross 
section related to u which is used in the Geant4 database is valid in an 
energy range from 1 keV to 10 GeV and taken from Seltzer and Berger 
(1985). The singly differential cross section related to 8,; i s based on the 
doubly differential cross section by (ITsai. 1974J : iTsai. 19771 ) and valid for 
very high energies, i.e., well above (1 — 10) MeV. But in the Geant4 database 
the dependence on the photon energy is neglected in this case so that it is 
actually a singly differential cross section describing Gj. 

Table [1] gives an overview of the available literature and data for total 
or singly, doubly or triply differential Bremsstrahlung cross sections; param- 
eterized angles or photon energies are given, as well as the different energy 
ranges of the incident electron. Furthermore, the table shows the atomic 
number Z investigated and includes some further remarks. 

For calculating the angularly resolved photon energy spectrum of TGF's, 
we need a doubly differential cross section resolving both energy and emis- 
sion angle of the photons; we need it in the energy range between 1 keV 
and 1 GeV for the small atomic numbers Z = 7 and 8. Therefore most of 
the literature reviewed here is not applicable. However, the Bethe-Heitler 
approximatio n is valid for atomic numbers Z < 29 and for electron energies 
above 1 keV (IShaffer et al.. 19961 ). How the range of validity depends on the 
atomic number Z is discussed in section 12.21 We therefore will use the triply 
differential cross section derived by Bethe and Heitler (1934) to determine 
the correct doubly differential cross. 
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Table 1: Available data for Bremsstrahlung cross sections. Besides the available information on total or singly, doubly or triply 
differential cross-sections, the range of validity of the incident electron energy and of the atomic number is given. If not stated 
otherwise, these are theoretical expressions. 
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1.4- Bremsstrahlung data used by other TGF researchers 

Carlson et al. (2009, 2010) use the Geant 4 database to simulate the pro- 
duction of Terrestrial Gamma-Ray Flashes. But Geant 4 does not supply an 
energy resolved angular distribution as it does not contain a doubly differ- 
ential cross section, parameterizing both energy and emission angle of the 
Bremsstrahlung photons (see Tabled]). We will compare the Bethe-Heitler 



approximation with Geant 4 in Appendix D 



Lehtinen has suggested a doubly differential cross section in his PhD the- 
sis (Lehtlinen, 2000) that is also used in a recently submitted manuscript 
by Celestin and Pasko. Lehtinen's ansatz is a heuristic approach based on 
factorization into two factors. The first factor is the singly differential cross 
section of Bethe and Heitler (1934) that resolves only electron and photon 
energies, but no angles. The second factor is due to Jackson (1975), it de- 
pends on the variable (1-/3^) [(1 — /3cos6i)^ + (cos6j — /?)] / (1 — /3cos6j)^, 
where /3 = \vi\/c measures the incident electron velocity on the relativistic 
scale. However, this factor is derived in Jackson (1975) only in the relativis- 
tic case, e.g., for incident electron energies above 500 keV. Furthermore, the 
calculation is classical and not quantum mechanical, and it is valid only if the 
photon energy is much smaller than the total energy of the incident electron. 
We will compare this ansatz with our results in Appendix E 



Dwyer (2007) chooses to use the triply differential cross section by Bethe 
and Heitler (1934), but with an additional for m factor parameterizing the 



structure of the nucleus (IKoch and Motz. 19591 ). This cross section depends 



on all three angles as shown in Fig. 1. If one is only interested in the angle 
Qi between incident electron and emitted Bremsstrahlung photon, the angles 
Gj and $ have to be integrated out — either numerically, or the analytical 
results derived in the present paper can be used. 

1.5. Organization of the paper 

In chapter [2] we introduce the triply differential cross section derived by 
Bethe and Heitler Then we integrate over the two angles G/ and $ to ob- 
tain the doubly differential cross section which gives a correlation between 
the energy of the emitted photon and its direction relative to the incident 
electron. Furthermore, we investigate the limit of very small or very large 
angles and of high photon energies; this also serves as a consistency check 
for the correct integration of the full expression. 
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In chapter in] we perform the same calculations for pair production, i.e., 
when an incident photon interacts with an atomic nucleus and creates a 
positron electron pair. As we explain, this process is actually related by some 
physical symmetry to Bremsstrahlung, therefore results can be transferred 
from Bremsstrahlung to pair production. We get a doubly differential cross 
section for energy and emission angle of the created positron relative to 
direction and energy of the incident photon. 

The physical interpretation and implications of our analytical results is 
discussed in chapter HI Energies and emission angles of the created photons 
and positrons are described in the particular case of scattering on nitrogen 
nuclei. For electron energies below 100 keV, the emission of Bremsstrahlung 
photons in different directions varies typically by not more than an order of 
magnitude, while for higher electron energies the photons are mainly emitted 
in forward direction. For this case, we derive an analytical approximation 
for the most likely emission angle of Bremsstrahlung photons and positrons 
for given particle energies. 

In chapter we will briefly summarize the results of our calculations. 



Details of our calculations can be found in [Appendix A| - [Appendix H 



Beyond that we provide a C++ code in Appendix I which can be used 
to create data tables of our analytical results for doubly differential cross- 
sections both for Bremsstrahlung photons and for pair production positrons 
for use in Monte Carlo codes. The C++ code can also be directly downloaded 
from the website of the journal. 
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2. Bremsstrahlung 

2.1. Definition of the process 

If an electron with momentum approaches the nucleus of an atom, 
it can change its direction due to Coulomb interaction with the nucleus; 
the electron acceleration creates a Bremsstrahlung photon with momentum 
k that can be emitted at an angle 0j relative to the initial direction of 
the electron. The new direction of the electron forms an angle / with the 
direction of the photon. The angle $ is the angle between the planes spanned 
by the vector pairs (pi,k) and (p/,k). This process is shown in figure 1. A 
virtual photon (allowed by Heisenberg's uncertainty principle) transfers a 
momentum q between the electron and the nucleus. Therefore both energy 
and momentum are conserved in the scattering process. 

The corresponding triply differential cross section was derived by Bethe 
and Heitler (1934): 
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Here Z is the atomic number of the nulceus, a fine ~ 1/137 is the fine struc- 
ture constant, h ^ 6.63-10"^^ Js is Planck's constant, h — h/27r and c 3-10^ 
m/s is the speed of light. The kinetic energy E^in^i/f of the electron in the 



kin,i/f 

initial and final state is related to its total energy and momentum as 
Ei/f = Ekin,i/f + rUeC^ = y^m^^^TpJ^ 



(2) 



where rUe ~ 9.1 • 10 kg is the electron mass. The conservation of energy 
implies 



Ef — Ei — hw 



(3) 
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which determines Ef asa function of Ei and Hoo. The directions of the emitted 
photon with energy hw and of the scattered electron are parameterized by 
the three angles (see Fig. 1) 

0. = <(p.,k), (4) 

0/ = <(P/,k), (5) 
$ = Angle between the planes (Pi,k) and (p/,k). (6) 

The differentials are 

dQi = sinBjfiBj, (7) 
dilf = sin e f dOf. (8) 

Furthermore one can get an expression for the absolute value of the virtual 
photon q with the help of the momenta, the photon energy hw and the angles 
dH) - ([6]). Its value is 



— q = — IpjI — Ip/I ~ j + 2|pi|— wcosBj — 2|p/|— wcosG/ 

+ 2|pj| |pj|(cos 6/ cos 6i + sin 9/ sin 9j cos $). (9) 

2.2. Validity of the cross sections of Bethe and Heitler 

The cross sections of Bethe an d Heitler ([1]) are valid if the Born approxi- 
mation ( iBethe and Heitler. 1934 ) holds 

Zc , . 

V » ^ (10) 

For nitrogen with Z = 7 and for oxygen with Z = 8, this holds for electron 
velocities |v^=7| ^ 15 ■ 10^ m/s and |v^=8| ^ 18 ■ 10^ m/s; this is equivalent 
to a kinetic energy of 

^2 



rriec' 2 ~- I 670 eV, Z 



= - » i 875 ev! Z = 8 " ^^^^ 

This means that incident electron energies above 1 keV can be treated with 
Eq. [H for lower energies, one cannot calculate with free electron waves any- 
more, but has to use Coulomb waves (Heitler, 1944; Greiner und Reinhardt, 
1995); in this case one cannot derive cross sections like ([T]) analytically any 
more. 
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2. 3. Integration over $ 

The easiest way is to integrate over the angle $ between the scattering 
planes first (see Fig. [T]) . For this purpose it is useful to redefine some 
quantities in the following way; therefore ([T]) can be written much more 
simply: 
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With (IT2]) - ([THD, Eq. ([T]) can be written as: 
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where Oj, i G {1, . . . , 4}, a and /3 still depend on 9/ and 9j. These integrals 
can be calculated with the help of the residue theorem which is reviewed 
briefly in Appendix A If i?(a;, y) : — M is a rational function without 
poles on the unit circle + y'^ = 1, then 
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The poles of the functions fi{z) in fl28p - 0311) are given by 
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It follows immediately that \zi\ > 1 and \z2\ < 1. For all residua one obtains 

Res(/i,^2) = -- J-— , (39) 

ResiKz,) = — ^=i=, (40) 

Res(/2,^3) = -, (41) 
m 

Res(/3,^2) = -^-^^^=7, (42) 



(v^ 



a 



2\3 ' 



Res(/4,Z2) = 040— ===— . (43) 



(v^ 



2^3 
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With the knowledge of these residua and using (1^ . the integral in (120]) can 
be calculated elementarily 



27rao 



duodVlidVl 



+ 



27r 



02/3 



as/? 



+ 



(44) 



2.4- Integration over Q f 

After having obtained an expression for the "triply" differential cross 
section, there is still the integration over By left. This calculation is mainly 
straight forward, but rather tedious. Using expression (jH]), it is 



_fa_ 

dudQj 



dQi 



2'Kan 



a 



+ 



27r 



— ai 



a 



+ 



|/32-a2| |/32-a2 



sin 9/. (45) 



Let's now consider the first integral of ( l45i) . If one inserts (fT2l) and (fT6l) . it 
becomes 



afc)/ sm 



E'i — cpi cos 9j 
27rAc2 



sin 



/ 



Ei — cpi cos Qi J Ef — cpfX 



dx- 



Ef - cpf cos 9/ 
1 



(46) 
(47) 



where the substitution x := cos 6/ was made in the second step. (H7|) is 
rather simple and yields 



d^yt sm 

a 



{Ei - cpi cos Qi)pf \Ef-pfC 



(48) 



Here "triply" really means the dependence on the photon frequency and two angles. 
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This was a quite simple calculation. All the other integrals can be calculated 
similarly, but with more effort. As another example let's consider the last 
integral. Before inserting f|T2|) . f lT3|) and f lTSj) one can define for simplicity 



dQ f — — sin 



(49) 



^1 := -Pj-P/-(-^l + 2-w|pi| cosBi, 
A2 := -2-w|p/| + 2|pj||p/| cosGj. (50) 



The expression /3 from Eq. ( JT3ll is then 

/3 = Ai + A2COS0/. (51) 

Thus the regularly appearing term — can be written as 

= (A^ + 4p^pJ sin^ Gi) cos^ 6/ + 2Ai A2 cos 9/ 

+ (Al-Apy^sm^Q,) (52) 

= cos^ 9/ + 2Ai A2 cos 6/ + Da (53) 

where the definitions 

Dl := A2 + 4p2pJsin2e„ (54) 
Dl := A^ - 4p,2pJ sin^ (55) 

have been introduced. 

By using ( IT2l) . ( JTSl) and ( |53l) . the last integral of (Hil) becomes 

27ra4a . „ _ 1671^^2^251^29.(^2^^2) 



^(/32 _ ^2)3 ^ - cpi cos Gi 

u 

X / d<df — , sm9 



□f cos ej + 2Ai A2 cos 6/ + niy{Ef - cpf cos 6/) 

(56) 

167rAp2p2 sin2 Q.^Ef + Ej) 



Ei — cpi cos 9j 

X / dx , (57) 

J v^(nfx2 + 2A1A2X + □i)3(Ej - cp^x) 



16 



where x = cos 9 / has been substituted again. 

This integration can be performed elementarily by finding the indefinite 
integral 



X (-DiE/a: + n2PfC+2AlAl(Ejx - c) + Ai A2 Da (E/ + P/ ca:) 



ln(^{Ef - pfCx)(DfEfX + D2PfC 



^(□fE2 + 2AiA2EfPfC + n-^pfc)^ 
AiA2(Ej! -\- pj!cx) 

^afEj + 2Ai A2-B/P/C + □|p_fc) ) , (S8) 



by inserting +1 and —1 as upper and lower limit, using ( IMI) and (1551) and 
simplifying. The integral in ( l56l) is then finally 



r 2naia 16n Ap^ p^f sin^ O i(E'f + E'i) 
/ de f — sin e f - ^ 



\/(/3^ - a^)'i ' Ei - cpi cos 0; 

2(A2P/c + AiB^) 



(-A| + A2 - 4p|p-^ sin2 ei)((A2-B/ + Aipjc)^ + 4m^c^p'^p^j, sin^ 0;) 



^((A2-B/ + Aipjc)2 + 4m2c<lp2p2 si„2 .)3 



X In 



^((E/ - cpi)Hp1p) sin^ ea-B/ " P/c) + (Ai - A2){{A2B/ + Aip^c) 
- ^Df + 2A1A2E fpfc + DIp/c)) ((iSj + cpf)(4plp) sin^ ei( + B/ -P/c) 
+ (Ai + A2)((A2B/ + Aipyc) - -ynfBl +2AiA2B/p/c + n|py.c))"^j . 

(59) 

All the other integrals can be calculated similarly where one always has to 
substitute x = cos 9/. With this technique the whole doubly differential 
cross section finally becomes 

(fa{Ei,u,Qi) A 
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with the following contributions: 



A| + ipjp) sin^ e, - ^A| + 4p2p2 sin2ei(Ai + Aa) + Ai Aa 
^A2 + 4p2p2 j.ijj2 0. I _a2 _ 4p?p2 j,;jj2 _ ^A2 + 4p?p2 ^^^2 ei(Ai - Aa) + Aj Aa 

cAa pfc^sin^Bi 2fi^iiJ^pyA2 



p^' (i?^ — cp^ cos 0^) (Si — cpi cos 0^)2 c(i!;i — cp^ cos 0^ ) (A| + 4p?p^ sin2 0^ ) 



(61) 



,a . _ l„|:gl±j^V (62) 

^3 = 



py(-Ei — cpicos0i) \^ E ^ — p j-c 
27rA 



^(Aa-B/ + Aipy:c)2 + 4,m^c1-p^pj sin2 Si 

In ^((-B/ + Pfc)(4,p^pj sin^ 0i(-B/ - p/c) + (Ai + Aa)((Aa-E/ + Aipyc) 

^(Aa-B/ + Aip/c)2 + 4m2c''pjp2 sin2 ©i))) ({B/ - p/c)(4pjp^ sin^ &i{-Ej - pfc) 

(Ai - A2){{AaB/ + Aipjc) - ^(AaB/ + Aipjc)2 + 4m2c4p2p2 ^j^a e^))) 

(A| +4p?p^ siii^ ©i)(Bj + Efpjc^) +pfc(2{A^ - 4p?pJ sin^ 0i)Bjpjc + Ai A2 (3£;| +Pjc2)) 
(A2E/ + Aipy:c)2 + 4m2c''p^p2 sin2 0; 

;(Aa-B/ + Aip/c) 



4Bfp|^(2(AaB/ + Aipyc)^ - 4m2c*pfp^ sin^ 0i)(AiEj + Aap/c) 



p y {Ej^ — Cpi COS 0^ ) 

{{A2EJ: + Aipjc)2 + 4m2c<lpfp^ sin2 0j) 
• + Ej) - 2n?u}'^p^ sin^ &iPfc{l 

{Ei - cpi cos©i)((Aa-E/ + Aipj:c)2 + Am'^c^p^p^^ sin2 ©;) 



Ip^m^c* sin^ ©iC-B? + -Bj) - 2fi^ij^p? sin^ 0ipjc(Aa-B/ + Aipjc) + 2R^i^^p j m^c^ ( Aa -B/ + Aipjc) 



(63) 

47rApj:c(Aa-B/ + Aipjc) 167rBf pjA(Aa-E/ + Aipjc)^ 

(Aa-Bj + Aip_f c)2 + 4m2c-*pjp^ sin2 ©; ((Aa-B/ + Aip/c)2 + 4m^c^p^p^^ sin2 G;)^ ' 



47rA 



(-A| + Af - 4p2p2 sin2 ©i)((Aa-B/ + Aip/c)2 + 4m2c''p2p2 sin2 0^) 
n^ui^pj 



Ei — cpi cos 0^ 

B/[2A|(A| - aI) + 8p?pJ sin^ 0i(A| + Af )] + py:c[2AiAa(A| - Af ) - 16AiAap?P/ sin^ ©i] 

A2 + 4p2p2 sin2 0^ 

2h^cj^ sin^ 0i(2Ai Aap/c + 2A|Ej: + 8pf sin^ ©i-B/) 
i?^ — cpj COS ©j 

2i;fp^{2(A| - Aj)(A2B/ + Aipj,c)^ +Spfp^ sin''0i[(Aj + A|)(i!;^ + p^c^) + 4Ai AaiSjpjc]} 



((AaE/ + Aip/c)2 + 4m2c'lpJpJ sin2 0j) 
p?pj sin^ 0i(B| + -E|)(Aap/c+ AiBj)" 



Ei — cpi cos ©^ 
lewEjp^ sin^ ©jA 



(B; - cpi cos 0i)2(-A| + Af - 4p2p2 sin2 ©;) ' 



(65) 
(66) 



Eq. f l60|) depends explicitly on Ei,u and Gj while Ef and p/ are functions 
of Ei and w through (|2]) and ([3]). ( 160|) is the final result of the integration 
of (II]) over $ and 0/ with the help of the residue theorem and some basic 
calculations. Now this result can be used both as input for Monte Carlo 
code and for discussing some basic properties of the behaviour of produced 
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Bremsstrahlung photons. 

Actually fl6Up is also valid for Bj = 0, as will be shown in the next section, 
but the simple way just to set Bj = in f l60|) will fail, especially for numerical 
purposes, because the logarithmic part in fl6T|) tends to ln(0/0) for Bj — )■ 
and so fails for numerical applications. Thus we need an additional expression 
for Bj = which has to be consistent with flBUl) . 

2.5. Special limits: Bj = 0, tt and Tiw — -Efcm,j 

For some special cases the integration of ([1]) over $ and B/ is easier. This 
information can also be used to verify ( 160|) by checking consistency and use 
them for Monte Carlo codes. 

2.5.1. Bi = or Bi = TT 

If one is only interested in forward and backward scattering, one can set 
Bj = or Bj = TT before integrating; then ^ becomes 



dojdVtidVtfd^^ (27r)2 |pj| u |q| 



X 



p^sin^Bf , „ „ „, 



{Ef — c\pf \ cosB/)2 



9 9 P f sin^ B f \ , X 

+ 2 h oj — - I (67) 

(Ef - c\pf\ cosB/)(^i =F c|pi| cosBi) J 

where the momentum q of the virtual photon can be written as 
— = — p^ — pj — ( —u j — 2 — IpjI cos B/ ± 2— cjIpjI ± 2|pj| |pj| cos B/. 

(68 

Here the upper sign corresponds to Bj = and the lower one to Bj = vr. 
As (167|) and (168|) do not depend on $ at all, the $ integration simply 
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gives a factor of 2tt, and fl^ becomes 



dujdVLidVL j 



27r 



P/ 



2 sin^ e 



/ 



+ 2 



(E;-c|p;|cOSe^)2 

pIsin^B/ 



2, ,2 



(i^^/ - c|p/| cos0/)(i?j =F c|pi| cos6i 



(69) 



Finally this expression has to be integrated over 0/ in order to obtain the 
doubly differential cross section. Similarly to the total integration of (H5|) it 
is convenient to define 



A f ^ \ 2 

Ai := - \pi =F -UJ 1 -pj, 

h 

A2 := -2-ujpf ± 2piPf 



(70) 
(71) 



where Ai,2 = Ai,2(ei = 0,7r),j G {1,2} , with definitions and ([50]). Eq. 
fl68ll can then be rewritten as 



and 

d^a 
dujdVti 



- q2 = Ai + AscosG/ 



(72) 



^'«}.ne^' Ip/I 1 



27r 



iPil w 



|p^|V sin^ 



(Ej - |p/|ccose/)2(Ai + A2 cosG/) 



4E2|p^|2gin2 0^ 



(E/- |p/|ccos0/)2(Ai + A2Cose/)^ 

-.2, ,2 



+ 



2nV|p/psin20j 



(Ei T c\v^\){Ef - |p/|ccose/)(Ai + AacosO/)^ 



sin 6i 



(73) 
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where the integration is rather elementary and can be performed by substi- 
tuting X = cos 6 / again. Thus ([72]) yields 



+ 



+ 



dudVti 

|p/hc- 



Z^a^h'' IpjI 1 




2|p/|c 



A2^/ + Ai|p/|c 



A2(A2^/ + Ai|p/|c 



2^,Ef\pf\c + ^2{Ej + Ip^I V) [Ej + |p;| 



(A2i?/ + Ai|p;|c)2 



Ef - |pj|c 

2, ,2 



{A^Ef + Ai|pj|c)2 (A2^/ + A,\pf\c){E, t c|p.|)A2 
8E!\pj\\A,Ef + A2IP/IC) (Ai - A2)(i?/ - \pf\c 



(A2S/ + Ai|py|c) 



(Ai + A2)(^/ + |p/|c 



+ 



+ 



2h'\pf\'u\2A^A,Ef + A2|p^|c + A2|p^|c) 



{A2Ef + A^\pf\c)\E,Tc\pi\)Al 



In 



In 



(A2^/ + Ai|p;|c)2(E, T c\p,\)c \Ef + \pf\c 



|P/|c 




(74) 



This expression is much simpler than ( l6Up . but only valid for 0j = or 
Qi = IT. Actually this expression can also be obtained by calculating the 
limit 0j — )■ or 0j — )■ TT in f l60|) : hence both expressions are consistent. The 



reader interested in details of this calculation, is referred to [Appendix B 



2.5.2. hw Ekin,i 

The other case which can be investigated easily is when almost all kinetic 
energy of the incident electron is transferred to the emitted photon, i.e.. 



E 



f 



Ei — hw = Eui„ i + rrieC — hw 



kin,i 2 

7- rripC 



(75) 



and 



Ip/I 




— mtc^ 



(76) 
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and consequently from Eq. 



- q 



kin.i 9 

> -Pi 



h 



-OJ 



+ 2—u\pi\ cos 0j 

c 



r ^ 5 . 



With these hmits it follows for the triply differential cross section ([T]) 



kin,! 



dudflidQ fd^ 



X 



(77) 
(78) 





Pf 


1 1 


Pi psin^ 9i 




|Pi 


a; 52 


XEi - c Pi cos9i)2 


{4E} + 6^ 


) + 


2 h^u 


2 pipsin^Gj 




{Ei - c Pi cos6i)i?/ 



(79) 



Actually (179|l depends neither on $, nor on 0j. Therefore 



27r 



/.*/si 



sinGjdBi 



47r 



(80) 



which leads to a very simple expression for the doubly differential cross sec- 
tion 



d'^a hLu^Ek 



dudili 





Pf 


1 1 




TT 


Pi 


u6^ 





IPipsin^ 6i 



{Ei - c|pi| cos6i 



:iAEj + 6c' 



+ 2 h'^oo'^ 





Pi 


^sin^e, 




- c\pi 


cos Qi)Ef 



(81) 



It seems that taking the limit hu — )■ -Efcin,i contradicts Eq. (fTTl) as the energy 
of the emitted electron should be larger than 1 keV ( 1111) . too. But ( 18T|) can 
be used for two purposes. 

As (ISTj) can be obtained, as well, by taking the limit \pf\ — )• in (1^ . 
the complicated expression f l60|) is verified. For further details the reader 



is referred to Appendix C Although this formula is strictly valid only for 
hw/Ekin,i — >■ 1 we will, however, see in section 11.1.41 that this cross section 
can also be used for calculating the most probable scattering angle for almost 
every photon energy. 
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3. Pair production 



Pairs of electrons and positrons can be produced if a plioton interacts 
with the nucleus of an atom. This process is related by some symmetry to 
the production of Bremsstrahlung photons. Bremsstrahlung occurs when an 
electron is affected by the nucleus of an atom, scattered and then emits a 
photon. So there are three real particles involved: incident electron, scat- 
tered electron and emitted photon. As the photon has no antiparticle one can 
change the time direction of the photon. For antimatter it is well known that 
antiparticles can be interpreted as the corresponding particles moving back 
in time. So one can substitute the incident electron by an positron moving 
forward in time. Thus by substituting emitted photon by incident photon 
and incident electron by emitted positron (due to time reversal and changing 
its charge) it is possible to describe pair production from Bremsstrahlung. 
Thus the emitted photon in the Bremsstrahlung process has to be substituted 
by the incident photon from the nucleus and the incident electron by the pro- 
duced positron. With these two replacenaents one gets the differential c ross 
section for pair production (IHeitler. 19441 : iGreiner and Reinhardt. 19 951 ) 



(27r)2/i 



IP+IIP-I 



sin^ e_ 



X 



|q| 



{E^ — c|p_| cos 



{AEl - cW) 



p?_ sin^ 9- 



{E^ — c|p+| cos 0_ 



2, ,2 



pi sin 



(4E! - cW) 



+ pi sin^ 6_ 



{E+ — c|p+| cos9+)(i?_ — c|p_| cos6- 



+ 2 





P+l 


P- 


sin 6+ sin 9_ 


cos $ 




C|P4 


_ COS 6+)(i?_ — c p_ 


cos 6_) 



{2El + 2El - c2q2) 



(82) 

where Z, a fine, h, h and c are the same parameters as in Eq. ([T]). a; is the 
frequency of the incident photon, E± and p± are the total energy and the 
momentum of the positron/electron with 



E. 



(83) 



Similarly to ([T]) there are three angles, 6± between the direction of the photon 
and the positron/electron direction, 6+ = <(p+,k),6_ = <(p_,k), and $ 
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is the angle between the scattering planes (p+, k) and (p_, k). The absolute 
value of the momentum of the virtual photon is 



Ef 










— > 




f ) 


Pi 


— > 


-Ph 




P/ 


— )■ 


p , 




u 


— )■ 


-w, 




e^ 


— > 


TT — 


e 




— ^ 








— ^ 


$ _ 


TT 



~q = ^ P+ ~ P- ~ + 2 p+ — w cos 9+ + 2 p_ — w cos 9_ 

\c J c c 

- 2|p+||p_|(cos6+cos9_ + sin9+sin9_cos<l>). (84) 
Algebraically one obtains (!82|) from ([1]) by replacing 

(85) 
(86) 
(87) 
(88) 
(89) 

+, (90) 

(91) 
(92) 

where the quantities on the left hand side are for Bremsstrahlung, and those 
on the right hand side for pair production. At the end one has to multiply 
with an additional factor to get the correct prefactor. With all the mentioned 
substitutions it is 

|P-(_| C 

Because of this symmetry the results for pair production follow easily from 
those for Bremsstrahlung. 

The direction of the positron relative to the incident photon is given by 
integrating (!82|) over $ and 9_. But this is the same exercise as to integrate 
([l]) over $ and 9/. Because of the symmetry between Bremsstrahlung and 
pair production one can take (|60l) and substitute (l85l) - (!92l) to obtain a 
doubly differential cross section 
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with the following contributions 



12 
13 



14. 
15 



Is 



(p)a(p) \ 



(A<P))2 + Ap\pt sin^ Si - y(A^''^)2 + 4p2^pi Sin2 e^A^*"' + A^f ' ) + A^^'a 



cA 



(p) 



2 2-2^ 
p^c sin fc?j 



p_ — cpj^ cos — cp_|_ cos ©i)-^ 

i) \1 



c(B+ - cp+ cosei){{A^''')2 + 4p2_p2 gijj2 e-) 



27rAc 



p_ — cp+ cos Gi) V 

27rA 



^{a')['>E_ + A<'''p_c)2 + 4m2c4p2_p2_ gi„2 Q. 



(S_ + P-c)C4p^pt sin^ ei(B_ - p_c) + (A^"* + A^''')((A^''^B_ + A^^'^.c) 



v(p)^ 



>(P), 



JP), 



- p_c) 



(A^P^iS. + a5''V_c)2 +4m2c4p2 p2 3i,j2 
(a'*'^ - A'f^XCA^P'iS. + a'^^p.c) - -^(A'^'iS. + A^f 'p_c)2 + 4m2c*p2^pi sin2 64))) 



P— — cp^ COS 0^ ) 

[((A<P')^ + 4p^p?. sin^ e,;)(-E?- + E_p_c) + p_c(2((A[''> f - 4p\p^_ sin^ ei)E_p_ 



. (p) •,2 , . 2 4 2 2 ■ 2 „ 1 
-A^ 'p_c) + 4m c p_j_p_ sin 



2 2 2 4.2 ^ /ct2 , rn2 \ fifc2 2 2 .2 ^ / a (p) 

" m c sin ©j(i?_(_ + — 2/i w p_)_ sin GiP_c(A2 



\(P). 



-c) 



[ - Sp^p- 

2ft^w^p_m^c^(A^''>B_ + A^P'p.c)] - cp+ cos ei)((A^''>B_ + A<'''p_c)^ +4m^c*p^p^ sin^ G^) 

iE\p1{2(i!.'-^'>E_ + A<P)p_c)2 - 4m2c*p2 p2^ sin^ e,){A<'''i3_ + a'^'p.c)" 



((A^P^B_ + A<'''p_c)2 + 4m2c4p2_p2_ ^^^2 Q.p 



47rAp_c(A^*''B_ H-A^P'p.c) 



167ri52p?.A(A(P'B_ +A<P'p_c)2 



(A^P'b. +Ai'''p_c)2+4m2c4p2^p2 sin^Si ((A, 

4-77 A 



(p), 



B_ + AfP'p_c)2 + 4m2c4p2^p^ sin2 Gi)2 



(-(a(''))2 + (A^-'j^ 



4pi 



sin2ei)((A<'''B 



JP), 



c)2 + 4m2c4p2^p^ sin2 Gj) 



h ijj p_ 



-[b_[2(a('')) = ((a(''))^ - (a('')) = ) +8p>i sin^G,((A(''))= + (a('") = )] 



-E_|_cp_|_ cos ©i 

P_c[2aWa<'')((a(''))^ - (A<''))=) - ISAf^'Af^'p^p^. sin^ G^]] [(A^") )^ + 4p5.p^_ sin^ G,]"' 



2R2^2 sin^ Gi(2A(f' A(f'p_c+ 2(A<''')2S_ + 8p^p?. sin^ Gii3_) 



— cp_|_ cos ©^ 



[2B>^_{2{{a(''')^-{a(''))=)(a('')j 



,(P).(P), 



_p_c]}][(A 



(P) , 



v(P), 



-\-AY'p_c) + 8p^p_ sm ©i[((A:^^0 + (^2 ) 

\2,^ 242 2 - 2^1—1 
i_c) + 4m c Sin ©^ 



Sp^pi Sin2 e^(El + s!.)(A('')p_c + aWe_) 



- CP_]_ COS G; 

leirEip^ sin^ GiA 



(B+ - cp+ cosGi)2(-(A<''')2 + (a(''^)2 - 4p2^pi sin2 G+) 
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with 



^^^HneC^ IP.IIP-I ^^0^^ 



(27r)2/i a; 
and with A^^'' , A2^^ defined as 

A;^^ := -|p+|2-|p_|2-0^^+2^^|p+|cose+, (102) 

h 

A^^ := 2-a;|pi| -2|p+||p_|cose+. (103) 
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4. Discussion 



4.1. Bremsstrahlung 

4-1.1. Angular distribution of Bremsstrahlung 

If electrons are scattered at nuclei, they can produce hard Bremsstrahlung 
photons with frequency u and direction Gj relative to the direction of the 
electrons. 

Figure |2] shows the doubly differential cross section fIBU]) for Bremsstrah- 
lung for several electron and photon energies. At first, the probability for 
generating photons decreases with increasing photon energy for fixed electron 
energy. This can be understood easily by applying ( IHT]) . As can be seen there, 
the doubly differential cross section grows linearly in the momentum of the 
electron in the final state which is equivalent to 

for high photon energies. So, if all kinetic energy is transferred from the 
electron onto the photon, the final momentum |pj| vanishes, and thus 

0. (105) 



dudQ 



For nonrelativistic electron and photon energies the scattering angle tends to 
be mainly equally distributed, i.e. the photons do not have a preference for a 
particular direction. When the photon energy increases, photons are mainly 
emitted in forward direction, but the ratio between forward and backward 
scattering is at least three orders of magnitude lower than for a relativistic 
electron. This case belongs to the classical case where the velocity is small 
compared to the speed of light. Namely, it is v/c\Eki„i=iokeY ~ 0.20 and 
non-relativistic equations will be enough to describe these phenomena. In 
the relativistic case {v/c\E^^^^,=iMeY ~ 0.94 and v/c\e^^^^,=ioomcV ~ 0.99999) 
the differential cross section becomes more and more anisotropic. Forward 
scattering is preferred to backward scattering although the maximal cross 
section does not lie precisely at Gj = as can be seen in figure |2] c). But 
the more the electron energy increases the more the maximum wanders to 
smaller angles, for example, it seems in figure [2] d) that the maximal emission 



is indeed for Gi=0. As mentioned in section and in [Appendix B[ formula 



( 16T|) cannot be evaluated directly at Gj = 0. However, for this purpose, we 



27 



>i 1e-46 

® 

3 



1e-47 
1e-48 



nco/E^i^i=99%"\^ 



20 40 60 80 100 120 140 160 180 
©;[■>] 



1e-49 
1e-50 
1e-51 



fico/E, 



kin,i' 



=50% 



20 40 60 80 100 120 140 160 180 



1e-44 

(<r~ 1 e-45 
E 

^ 1e-46 
1e-47 
1e-48 
1e-49 
1e-50 
1e-51 
1e-52 
1e-53 



m 
o 



= 10 keV 






fico/E|(|p 1=1 % 




fl<0/E|,|n|=10% 




flO)/E|(ini=56% 




fi<o/E,i„ 1=99% ; 



c) 



kin.i 



20 40 60 80 100 120 140 160 180 

1 MeV 



1e-40 
1e-42 
1e-44 
1e-46 
1e-48 
1e-50 
1e-52 
1e-54 
1e-56 
1e-58 
1e-60 



150 keV 



E 



® 
S 



LU 



Ii<o/E|^iP i=1% 



fio)/E|(i„i=10% 

Fico/Eki„ 1=50%. 



fio)/E|(ln 1=99% 



20 40 60 



80 100 120 140 160 180 



kin,i 



100 MeV 



Figure 2: The doubly differential cross section (P<T/{dujdil.i){Ekin,i,^^,'di) for 
Bremsstrahlung {Z = 7) versus the scattering angle 0j between emitted photon and 
incident electron. The electron energies are a) Ekins = 10 keV, b) Ekin,i = 150 keV, c) 
Ekin,i ~ 1 MeV and d) Ekin,i = 100 McV. In each plot the the photon energy tuj amounts 
to 1%, 10%, 50% and 95% of the kinetic energy of the incident electron. 
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Figure 3: The doubly differential cross section (P a j (diod^i) as in figure [2] for a smaller 
angular range 
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Figure 4: The ratio between the doubly differential cross section for backward scattering 
(Oi = 180°) (P(j{Ekin,i,uj,&i = 180°) / {dudfti) and the maximum of this cross section 
max (c?'^(T(£'fci„^i,a;, 0i)/(c?aj(ii7i)) vs. the kinetic energy of the incident electron for differ- 
ent ratios between photon and electron energies in a a) linear and b) logarithmic scale for 
Z = 7 . 
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Ze Ze 




Figure 5: a) In the rest frame of the electron where the nucleus is moving instead radiation 
is emitted isotropically with a small-angle deflection, b) If, however, one transforms the 
situation into the rest frame of the nucleus where the electron is moving, most of the 
radiation is emitted in forward direction relative to the direction of the electron. 



derived (17il) which is vahd for 9j = and 9j = vr. Figure [3] shows again (160!) 
for two relativistic electron and different photon energies but for a smaller 
range of angles. It shows in more detail that the angle of maximal scattering 
is small, but not 0. 

Figure H] shows the ratio between the cross section for backward scat- 
tering and forward scattering. It can be clearly seen that the tendency for 
backward scattering decreases for increasing electron energy. The lower the 
electron energy becomes, the more forward and backward scattering become 
similar and in general, the scattering tends to be isotropic. Only for ratios 
between photon energies and electron energies close to 1, forward scattering 
is preferred for the whole range of energies, but still decreases with increasing 
electron energies. 

4- 1-2. Relativistic transformation 

The tendency of forward scattering in the case of relativistic incident 
electrons can be understood by applying the laws of relativistic transforma- 
tions. Imagine a non-quantum field theoretical description of Bremsstrahlung 



(I Jackson. 19751 ). If one regards an inertial system in which the incident parti- 
cle is at rest (Fig. [5]a) ), radiation is emitted isotropically with a small-angle 
deflection. If the physical laws for this process are relativistically transformed 
into the laboratory system where the nucleus is at rest and the electron mov- 
ing, most of the radiation is emitted in forward direction relative to the 
electron direction (Fig. [5]b) ). Because this transformation is valid for a 
non-quantum field theoretical, relativistic electron, it must also be true for a 
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Figure 6: The doubly differential cross section (Pa/{dujdili){Ei:ini,uj,Qi) for 
Bremsstrahlung (Z = 7) versus the scattering angle Qi between emitted photon and 
incident electron. The electron energies are a) Ekin,i — 400 keV, b) Ekin,i ~ 511 keV. 
In each plot the the photon energy huj amounts to 1%, 10%, 50% and 95% of the kinetic 
energy of the incident electron. 



relativistic quantum theoretical description, therefore we see that the forward 
scattering of photons can simply be explained as a result of the relativistic 
transformation. 

The forward scattering can moreover be understood by using the conser- 
vation laws of energy and momentum. They predict that photons have to be 
scattered in forward direction if electron and photon energy are high. The 



interested reader is referred to Appendix G 



Although figure [2] shows that the maxima of the doubly differential cross 
section form with increasing electron energy, it is difficult to determine in 
these plots when these maxima really start to be generated clearly. Figure 
in] shows the doubly differential cross section in dependence of the incident 
electron energy for -Efcj„_j = 400 keV and -Efcm.i ~ 511 keV when the kinetic 
energy is equal to the rest energy. For 400 keV and for hw / Ekin,i = 0.01 the 
cross section for forward scattering is already two orders of magnitude larger 
than for backward scattering, but a clear maximum cannot be seen. How- 
ever, for the same kinetic energy, but for tvjj/Ekin,i = 0.95 there is already 
a clear maximum formed. But if the kinetic energy grows up to 511 keV 
which is equal to the rest energy of the electron, there is even a maximum 
for hw = 0.01Ekin,i- This can be expected due to the relativistic transforma- 
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tion. If -E'fcm,j <^ fneC^, then the photon emission is relatively isotropic. But 
if the kinetic energy is approximately equal to the rest energy of the electron, 
relativistic laws are valid. Especially for Ekin,i = rrieC^ 

- = ^^ 87%, (106) 
c 2 

therefore the electron has to be treated relativistically and clear maxima 
close to Gj = form for every possible photon energy. 

4-1.3. Dependence on the energy of the emitted photon 

Figure [2] shows that for both slow and relativistic electrons the doubly 
differentical cross section also varies with the photon energy for fixed electron 
energies. For fixed electron energy, low er photon energies are more likely. 
Moreover, the probability of photon generation is larger for low energetic 
electrons, but this is only true on average. Figure [7] shows the doubly differ- 
ential cross section in another way. Now the photon energy is fixed and the 
electron energy differs within one plot. For all cases it is more likely that low 
energetic electrons create photons than relativistic electrons do. This can be 
understood if one goes one step back to the triply differential cross section 
(HI). Considering that due to ([2]) and (E]) |pi| ~ Ei and jqp ~ jp^p ~ Ef, (H]) 
behaves like 

dujdVL.dVLfd^ ^ Ef' 

therefore ([T]) decreases if the energy of the incident electron increases. And 
because of the fact that integration preserves monotony, the doubly differ- 
ential cross section ( 160|) has to decrease for increasing electron energy as 
well. But for small angles, i.e. for forward emission of photons, the proba- 
bility rapidly increases for relativistic electrons and exceeds the probability 
at small electron energies. 

4.1.4- The most probable scattering angle 

Figure [2] also shows that the angle for which maximal scattering takes 
place, is rather independent of the photon energy. Hence, one can use f lHTj) 
to determine a formula for that scattering angle. Actually this derivation 
leads to a quartic equation which can, however, be approximated for small 
angles, i.e. 0j < 20°, through a quadratic equation. The reader is referred 
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Figure 7: The doubly differential cross section d'^a-/{dujdni){Ekin^i,uj,'di) for 
Bremsstrahlung {Z = 7) vs. the scattering angle 6j for several electron and photon 
energies. In each panel the photon energy tuv is fixed and the cross section is plotted for 
various kinetic energies of the incident electron. 
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Figure 8: Qi for maximal scattering vs. incident electron energy in a a) semilog and b) 
loglog plot for Z = 7 . Besides (|108p for huj — 0.9999Ekin,i, the exact solution of the 
quartic equation and various data for different huj/Ekin,i are shown. 



to Appendix H for the detailed calculation. The solution of the quadratic 
equation reads 



with 



1 |pi 



-^(4E| + 5oc2)-^(i?.-c|p,|) 



AEj + 5oc2 + - c|p,|) - |p,|5oC - f c|p,|5, 



(108) 



' -co] + 2-uj\pi 



c 



huj Ekin,i, e.g. hw = 0.9999Ekin,i and 



c 



(109) 



Pi 



kin,i 



kinA 



+ 2m, 



;iio) 



Figure M shows f llOSp and manually extracted values for 0, for different pho- 
ton energies. It shows much better than figure[2]that 0j is rather independent 
of the photon energy for relativistic electron energies. Besides fllOSp . the so- 
lution of the quartic equation, is shown. Moreover, it shows that (11081) gives 
a good approximation for those angles Gj for which scattering is maximal. 
Actually, we see that the exact solution describes the angle for maximal scat- 
tering better, especially for low energies, but for high energies both curves 
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fit very well. 

By inserting Ekin,i = huj/0.9999 into (I1U8I) one obtains a formula which 
relates the photon energy to the most probable scattering angle. 

4-2. Pair production 

4-2.1. Basic properties of pair production 

We now proceed from Bremsstrahlung to pair production. One photon 
with energy hw creates two particles, namely an electron and a positron, 
both with rest energy rrieC?- Therefore 



follows for the kinetic energies of these two particles. Thus the photon en- 
ergy has to be /zo; > 2meC^ 1.022 MeV for pair production and the kinetic 
energy of the particles is bounded as Ekin,± < JtjjJ — 2me(?- Figure [9] shows 
the doubly differential cross section flMj) for different photon and positron 
energies. Forward scattering is preferred, but there is almost no case now of 
more isotropic scattering . This results from the fact that the incident parti- 
cle is a photon and these have always to be treated relativistically. For very 
highly energetic photons, e.g. 50 MeV and 100 MeV, there are clear maxima 
for forward scattering. Additionally positrons with higher energies are less 
probable than positrons with lower energies which is the same as for photons 
emitted due to Bremsstrahlung. For energies /to; < 50 MeV , however, the 
maxima are > 5° , but forward scattering is still preferred. However, there is 
a new phenomenon which has not occurred for Bremsstrahlung or for highly 
energetic pair production. The production of Bremsstrahlung photons with 
higher energies has always been more probable than for small energies. But 
figure [9] also shows that the probability for highly energetic positrons can 
exceed the probability for lowly energetic positrons, at least for small an- 
gles. Figure [TD] shows flMl) for fixed positron and different photon energies. 
Again positrons which are generated with high velocities prefer to be scat- 
tered forward while this tendency vanishes if the positron energy is very low. 
This can be traced back to the relativistic behaviour again. If a positron is 
very energetic, it has to be treated relativistically and the relativistic trans- 
formation leads to forward scattering (this is the same explanation as for 
Bremsstrahlung). We also see that the creation of positrons is more likely 
for highly energetic photons. 




hjjj — 2me(? 
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Figure 9: Doubly differential cross section (Pa{E_^.,uj, Q^)/{dE^dQ^) vs. 6+ as an angle 
between incident photon and created positron for Z = 7: The cross section is shown for 
a.) hoj = 5 McY, b) fiw = 10 McV, c) hu = 50 McV and d) hu = 100 MeY. In each 
panel different positron energies relative to the available photon energy huj — 2meC^ 
are plotted. 
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Figure 10: Doubly differential cross section (Pa{E-^-,bJ,Q-^-)/{dE+dil+) vs. 0+ as an 
angle between incident photon and created positron for Z ~ 7: The cross section is shown 
for a) Ekin,+ = 150 keV and b) Ekin,+ — 1 MeV. In each panel the photon energies are 
huj = 5 MeV, 10 MeV, 50 MeV and 100 MeV. 



4-2.2. The most probable scattering angle 

As for Bremsstrahlung one can get a simple formula for the preferred 
direction. Performing the same calculation as for Bremsstrahlung one obtains 



with 



and 



6, 



(p) 
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|P^ 



nio 
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(E+-c|p+|)) 



E_ 



c|Ph 



5^P'^ := -I 



'112) 
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- I -cu 



2 p+ — 
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hu — 2meC^ 

Figure ITT] shows that f lll2p is a good approximation for B+ for high photon 
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Figure 11: 0+ for maximal scattering vs. incident photon energy in a a) semilog and b) 
loglog plot for Z = 7 . Besides (|112p for Ekm^+/{fiaj — 2raf.c?) = 0.9999 various data for 
different Ekin^+/{fiui — 2m eC^) are shown. 



energies and high ratios between photon and positron energy. The smaller 
the ratio between photon and positron energy, however, is, the worse flll2p 
becomes for low photon energies. If the photon energy is larger than 50 MeV, 
relativistic positrons are created; therefore forward scattering takes place and 
9+ can be calculated with (I112p . 
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5. Conclusion 



We have reviewed literature relevant for Bremsstrahlung in Terrestrial 
gamma-ray flashes (TGFs) (Bethe and Heitler, 1934; Heitler, 1944; Elwert 
and Haug, 1969; Seltzer and Berger, 1985; Shaffer et al., 1996; Agostinelle 
et al., 2003). Focussing on atomic numbers Z = 7 (nitrogen) and Z = 8 
(oxygen) and an energy range of 1 keV to 1 GeV, no good parametrization 
of an energy resolved angular distribution in the form of doubly differential 
cross section is available. The theory of Bethe and Heitler covers this energy 
range for Z = 7, 8, but it parameterizes the direction of the scattered elec- 
tron as well; therefore we integrated their triply differential cross section to 
obtain the correct energy resolved angular distribution for Bremsstrahlung 
and pair production. Other authors (Lehtinen, 2000; Dwyer, 2007; Carlson 
et al., 2009, 2010) used different approaches, as discussed in the introduction. 
They use singly or triply differential cross sections which do not give a direct 
relation between the photon energy and the direction of the photon relative 
to the r notion of the electro n. As positrons are created within a thundercloud 
as well (iBriggs et al.. 20101 ). we used a symmetry between the production of 
Bremsstrahlung and the creation of an electron-positron pair both in the field 
of a nucleus to obtain a cross section which relates the energy of the created 
positron with its direction. 

We have seen that emitted Bremsstrahlung photons are mainly released 
in forward direction if the electron which interacts with the nucleus has such 
a high energy that it has to be treated relativistically. For lower energies 
scattering tends to be more isotropic. For the case that almost all kinetic 
energy of the incident electron is transformed into photon energy, we derived 
an approximation for the most probable photon emission angle as a function 
of the incident electron energy and of the photon energy. The expression 
is valid for all ratios of photon over electron energy if the electron motion 
is relativistic. So, when photons have been created within a thundercloud 
or discharge, they are mainly scattered in forward direction as long as the 
electrons move relativistically, i.e. if their kinetic energy is at least as large 
as their rest energy. 

Similar results hold for pair production. Next to the doubly differen- 
tial cross section we derived a simple approximation for the most probable 
positron emission angle for the case that the photon energy is larger than 
10 MeV (for ratios between the kinetic energy of the positron and available 
photon energy down to 25%) or than 100 MeV (for ratios lower than 25%). 
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We have seen that for very highly energetic photons that almost all positrons 
are scattered in forward direction. If, however, the photon energy decreases, 
the probability of forward scattering decreases as well. Instead the maxi- 
mal cross section can be found at G+ ~ 90° for low ratios between £"+ and 
huj — 2meC^ and is, beyond that, symmetric to this angle. 

Our analytical results for the doubly differential cross-sections for Brems- 
strahlung and pair production are also supplied in the form of two functions 
written in C++. In this form the functions can be implemented into Monte 
Carlo codes simulating energetic processes like the production of gamma-rays 
or electron positron pairs in thunderstorms. 
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Appendix A. The residual theorem to calculate integrals with trigono- 
metric functions 



In this appendix the method how to calculate integrals of the form 



2-K 



/2(cos$,sin$)(i$ 



(A.l) 



shall be discussed where R{x,y) : \ {x,y G M|?/ = — x"^} — )■ M is a 
rational function without poles on the unit circle + = 1. But before 
explaining this method let's briefly review some general facts about residua. 

Appendix A.l. The residual theorem 

Let / : C D / — )■ C, 2 H- be a holomorphic function and T : [a, b] — )■ 

C, if: r(t), a closed curve in the complex plane. Then one can calculate 
complex curve integrals via 



J fiz)dz = 2niY^Res{f,Zj] 



(A.2) 



where the sum has to be taken over all poles Zj of / and complex curve 
integrals are deflned as 



I f{z)dz := I fim) . ^{t)dt. 



(A.3) 



The residuum of a pole zj can be calculated via 



Res{ f,Zj] 



lim 



d 



n-l 



{n — 1)\ z^zj dz"' ^ 
where n denotes the order of the pole. 



{z - Zj)f{z) 



(A.4) 



Appendix A.2. Integral with trigonometric functions 

With the help of the (IA.2P one can simply perform the integration of 
(lA.ip . For that purpose deflne 



tz V 2 V z 



1 

2i 



z — 



(A.5) 
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and choose the unit circle 

r(t) = e^*,t G [0,27r] (A.6) 
as closed curve; hence ( IA.3P becomes 



f{z)dz 



I {I + ' h - ^~')) (^-'^ 



2n 

= j {d' + e-*) , 1 (e^* - e-*) j t/t (A.8) 



= j R{cost,smt)dt (A.9) 



where the identities cost = 1/2 (e** + e~**) and sint = l/(2z) (e** — e~**) were 
used in the last step. 

Finally with ( ]A.2p and ( lA.Qp one gets a simple formula to calculate ( lA.ip : 



r 

2n 



2n 



/ i?(cos $, sin <^)d<^ = ^ Res(/, z) (A. 10) 



with / being defined in f lA.Sp . 



Appendix B. The doubly differential cross section for 0j = and 

In order to get ( 17^ from ( l60l) it is rather straight forward to set 0^ = 
or 6j = vr. Especially it is 

Ai(e, = 0,7r) = Ai, (B.l) 
A2(ei = 0,7r) = A2. (B.2) 

But there is one case which should be considered a bit more thoroughly. 
This regards the logarithm in ( 1611) . For 0j = tt it is A2(0 = tt) = A2 = 
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—2pf{h/c u + Pi) < 0; thus IA2I = — A2 and 



-A|-4p2p2 gin2 ei-^A|+4p2p2 sin2 e,(Ai-A2)+AiA2 



IA2 



In 



— ^In 
A2 \ Ai - A2 



A2 + 4p,2p/ sin^ 



Al- |A2|(Ai + A2) + A1A2 
-A2' - |A2|(Ai - A2) + A1A2, 
Ai + A2' 



(B.3) 

(B.4) 
(B.5) 



which is a very simple calculation. However, for 0j = it is A2(0i = 0) = 
A2 = —2pf{h/c u — pi) which can be both negative or positive depending on 
values of Pi and /i/c a; If A2 < then equations fIB.Sp - fIB.Sp are valid again. 
If A2 > 0, however, it follows for the argument of the logarithm 



Ai + 4p2p2 sin2 Q. _ /^2 ^ 4p2p2 gj,^2 q^(^^ _^ _^ ^^^^ 



-Ai - 4p2p2 gin2 e. _ + 4p2p2 sin2 Q.(^^^ _ ^ ^^^^ 

Aj- A2(Ai + A2) + AiA2 _ 
-Ai - A2(Ai - A2) + A1A2 ~ ■ 

Hence it is necessary to use the rule of L 'Hospital: 



(B.6) 

ei=o 

(B.7) 



lim ■ 



lim 



Ai + Apfp) sin^ Qi - J^l + Apjp) sin^ ei(Ai + A2) + A1A2 



-Ai - 4p2p2 sin2 - ^A2 + 4p2p2 gj,^2 Q.^^^ _ ^ ^^^^ 

(B.8) 



In 



-8p2p2 sin e, cos e, - ■ ^Pip) sin 6,, cos 6 

Ai- A2 

Ai + A2 



(B.9) 
(B.IO) 
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With f lRTOj) the whole hmit yields 



. A2+4p2p2 si„2 0_^^2^_4p2p2 si„2 ei(Ai+A2)+Ai Aa 

-A2-4p2p2 gin2 ei-^A|+4pfp2 sin2 e^(Ai-A2)+AiA2 



e»=o 





(B.ll) 



(B.12) 



which is and has to be identical with (IB.Sp . So in both cases, A2 > and 
A2 < 0, ( 1B.5|B.12p are generated by setting 9j = 0, vr; therefore one does 
not have to distinguish between the these cases in (17^ . 

But it is of importance to mention that due to (IB. 711 one can get numerical 
problems if one only implements ( l60l) and wants to calculate the doubly 
differential cross section for 9j = 0. Thus it is useful to distinguish for 
6i 7^ and 0j = and to use f l7i|) instead for the latter case. 

For the rest of limiting forward and/or backward scattering it is, however, 
straight forward to insert 6^ = 0, vr and thus can deduce from (I6OI) with 
the additional help of fIB.lOp . 



Appendix C. The doubly differential cross section for huj — >• E^in^i 

There are three contributions from fl60|) which lead to flHTl) in the limit 
hu Ekin,i |p/| 0: 



(Ei - cpi cos ei)^(Af - - 4pfpj sin2 G;) 



2irAp?c^ sin^ G^ 



(Ei - cp, cos G,)2 ^A| + 4p2p2 ^i^2 



A| +4pfpJ sin^G; - ^Al + 4p?p2 sin2 Gi(Ai + Aj) +A1A2 
-A2 - 4p?p2 sin2 Gi - ^A2 + 4p2p2 siii2 Gi(Ai - A2) + AjAa 



CC.I) 



(C.2) 



4,Trh''u}''pf sin-' G; A 
Ei — cpi cos @i 



2Ai A2PJC + 2A|Bj, + 8p^p^ sin2 G^Ej: 



(-A| + Af - 4p2p2 sin2 Gi)((A2-B/ + Aipjc)2 + 4m2c*p?pJ sin2 G;) 



Pfc{A2Ej + Aipj:e) 



^((A2Ef + Aipj,c)2 + 4m2c<lp2p2 sin2 Gi)3 

In |^({-E/ + Pfc)(4p^pj sin^ Gi(B/ - pj:c) + (Ai + A2)((A2-B/ + Aipjc) 
^J{A2Ef + Aipjc)2 + 4m2c'4p2p2 sin2 Gi))) ((B/ - p/c)(4p^p^ sin^ Bii-Ej - pfc) 
(Ai - A2)((A2B/ + Aipjc) - ^{A2Ef + Aipjc)2 + 4m2c''p2p2 sin2 Gi))) 



(C.3) 
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while all other integrals which appear in f l6Up cancel each other (which will 
be shown in an example later). It can be verified easily that 



lim Ai 
lim A2 

Pf^O 



5, 



with A2 ~ |p/| 



(C.4) 
(C.5) 



according to definitions fj49l) . fl50l) and fl77|) . With these limits the behavior 
of ii for small pj can be calculated in a straight forward way: 



lim Li 



IGirAEjp^ sin^ 0^ 



(C.6) 



Pf^o [Ei - cpi cos ©4)252 
For ( 1C.2I1 and (1C.3[) . however, there is more effort to be invested. As it is 

sin^ 9, 



Al + Apfpj sin^ and In ^ (^A^ + Apjp) sin^ 6^ - ^ Al + Apfp^j i 

X (Ai + A2) + AiA2)/( - A2 - Apfpjsin^Qi - ^Al + Apfpjsin^ Qi{Ai 

— A2) + AiA2j^ — > for pf — 7- 0, one has to use the rule of L'Hospital. If 
one rewrites 

A2 = ^Pf (C.7) 



with 



this rule leads to 



lim 

p/->0 



^ := -2-u + 2piCosQi 
c 



JA^ +Apfpj sin^e 



X In 



Al + Apjp) sin^ Gi - J Al + Apjp) sin^ ei(Ai + A2) + A1A2 



2 

7' 



Al - Apjp) sin^ Qi - J Al + Apfp) sin^ ei(Ai - A2) + Ai A2 / 



(C.9) 
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thus 



_ AnAp'^c^ sin^ 6^ 
P/^o ''^ {Ei - cpi cos 9i)25 ■ 

The hmit of ( \C.3\i can also be calculated by using f lC.7p . It is 



(C.IO) 



4TTh-'ui''p^ sin^ BiA 
Ei — cpi cos 



■e^Ef) 



pJ(-A| + - 4p2 sin^ Bi)((9Ei + Aic)2 + 4>7 



pjc(9Ej + Aic) 



pJ((*Ej: + Aic)2 + 4m2c''p| sin^ 0;) 
1 



Pf^(<iEf + Aic)2 + 4m2c-lp2 gin2 0.) 
X In ^((Bj +pj.c)(4p2p2 j,i„2 Q^jg^ _ p^^j _^ ^ + Aic) 

- pf^{-^Ei + Aic)2 + 4m2c4p2 gin2 ej)^ (^(^e j - p j: c) (4p^p^. sin^ 0^ (-Bj -P/c) 
+ (Ai - *p|0(P/(*E/ + Aic) - pf^{mEj + Aic)2 + 4m2c'lp2 sin2 0;)))"^ j 



(C.ll) 



While can simply be reduced in the fractions, one has to use the rule of 
L 'Hospital again for the logarithmic part because it is 
PJ^J{^!Ef + Aic)2 + Arn^d^pl sin^ 9j and the logarithm — )■ for — )■ 0. 
Its limit is 



X In |^({-B/ +Pfc)(ip^pj sin^ 0i(-B/ " P/c) + (Ai + 9pf){pf(-i!Ef + Ajc) 

- Pf^i'iEf + Aic)2 + 4m2c''p2 sin2 ({-E/ - P/c)(4p;pJ sin^ @ i(- E f -p/c) 

+ (Ai - *pj)(pj:(*B/ + Aic) - pj-^{mEi + Aic)2 + im'^c'ipf sin^ 0;))) 



(C.12) 



thus the whole limit yields after some further calculations 

Sirffoo^p^sin^QiA 
lim ^3 = — - — . 

P/^o [Ei - cpi cos QijO^Ef 



(C.13) 



Finally, if one inserts (JH]), the sum of (JCJ]), f laTOj) and f lal3|l leads to (|8T]). 

All other terms which appear in fl60|) vanish. For this purpose one should 
regroup all terms according to their origin. As an example let's consider the 

three contributions which have arisen from f f °^ '^°^'^„ d^dQ r. For this 

0j=O <I>=0 
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integral it follows 



/ 7- 

Gf =0 * = 



ct2 cos ^ 

c 

COS ^ + /3 



{Ei — cpi cos &i)pf 
2tx Ac^ 



Ei — cpi cos j 



^A| + ip^p'j sin2 &iPfC 



A| + 4p?p^ sin^ ei - ^A2 + 4p2p2 gi„2 ei(Ai + A2) + A1A2 
-A| - 4p?p2 sin2 e; - ^A2 + 4pfp2 sin2 G^Ai - A2) + AjAj 
A2-E/ + Aip/c 
PfC^(A2Ef + Aipfc)'^ + 4m2c*pjp2. sin2 q2 

In ^((-B/ + P/c)(4p^p^ sin^ @i{Ej - p/c) + (Ai + A2)((A2-B/ + Aipjc) 

^{A2Ef + Aip fCp + 4m2c4p2p2 sin^ 0^)) ((B/ - P/c)(4p^p^ sin^ SiC-Bj - pjc) 

(Ai - A2)((A2Bj + Aip^c) - yl{A2Ef + Aip/c)2 + im'^c^p^pj ^in^ e^))) j 

2c / 2\ , , 2 1 

__ + *(__)+(*B,+.c) _J 



- cpi cos 0j 



where we have used (1C.9I) . ( 1C.12I) and 

hmlln (^l±Pi^ 

Pf^opf \Ef-pfC 



2c 



(C.14) 



(C.15) 



in the limiting step. Of course, this term has to vanish because 02 ~ Pf, but 
the concrete calculation after having integrated over $ and 0/ is much more 
complicated. Therefore we have just given an example here. Similarly, all 
other terms cancel so that only the limits of Li, i G {1, 2, 3}, stay. 



Appendix D. Comparison with Geant 4 

As mentioned in the introduction cross sections for Bremsstrahlung can 
be found in the Geant 4 database (Agostinelli et al., 2003).(j(i?A;m,j). Geant 4 
contains data for the total cross section a, the singly differential cross section 
dajdbj and a singly differential cross section da/dVLi depending on 0j, but 
not on uj. Figure [D.12I shows the comparison of the singly differential cross 
section da/du{Ekin,i,i^)', here the singly differential cross section for Brems- 
strahlung which can be obtained by integrating over $, Bj and G/, reads 
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Figure D.12: Singly differential cross section da/doj (on linear scale) for Z — 1 (N2) from 
the between Bethe - Heitler approximation (|D.1|) and the data of Geant 4 as a function 
of incident electron energy Ekin,i and photon frequency w. In the different panels, the 
energy of the emitted photon is fixed to tiuj = 1 keV (a), 150 keV (b), 1 MeV (c) and 10 
MeV (d). 
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fiBethe and Heitler. 1934 : 



Heitler. 19441 ) 



da 
duj 



+ L 




2E,E 



+ 



+ 



^3p3 c^pipj- 



3 c^PiPf c^p^p j 



mlc^hw / EiEf + cV - EiEf + cVf 2hujEiEf 
+ TT^^ ^^r^5 ^-^5 + 



2c'^PiPf 



(D.l) 



with 



EiEf + PiPfC^ — mlc* 

rriec'^fioj 
Ej + piC 



2 log 
2 log 
2 log 

137 



(D.2) 
(D.3) 
(D.4) 
(D.5) 



where tq k, 2.80 ■ 10~^^cm is the classical electron radius. 

The Geant 4 database uses theoretical formula for both the low and high 
energy range which also includes the shielding du e to the outer electrons 
which is actually not included in ([1]). ( ID.ip and (lAgostinelli et al.. 20031 ) 
vary from each other for all related electron and photon energies. Fig- 
ur e ID. 121 shows a qualita tvely different behaviour as for the based data 



( ISeltzer and Berger. 19851 ) mainly contains data for high atomic numbers 
which are not appropriate to describe Bremsstrahlung in air. Schaffer et al. 
(1996), however, argue that the theory of Bethe and Heitler is good for an 
energy range of 1 keV - 1 GeV as long as Z < 29, hence appropriate for the 
molecules of air. 



Appendix E. Comparison with Lehtinen (2000) 

Fi gure lE . 1 3 1 shows the comparison of f l60|) and the doubly differential cross 
section used by Lehtinen. In all plots the photon energy liw amounts 50% 
of the incident electron energy. Qualitatively there is a good agreement, but 
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Figure E.13: Comparison of the product ansatz from Lehtincn (2000) with our result (jBO]) 
of the integration of ([T]) for different electron energies (Z — 7): doubly differential cross 
section versus the scattering angle 8^ between incident electron and emitted photon. The 
ratio between the photon energy huj and the kinetic electron energy Ekin.i is fixed to 50%. 
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Figure F.14: Contribution of ([6T|) - ([66|) to (|60|) in a semilog plot for different electron and 
photon energies {Z — 7). 



both curves deviate some orders of magnitude. According to (IJackson. 19751 ) 
Lehtinen uses a product ansatz for the angular and the freq uency part; her e 
the angular part is a non-quantum mechanical expression (IJackson. 19751 ). 
The product ansatz is said to be valid for very high electron energies and 
very small photon energies; therefore this cross section does not cover the 
whole energy range needed to simulate Bremsstrahlung in air. 



Appendix F. Contribution of the integrals 

As equation f l60|) is rather complicated, it is interesting to see which terms 
have the most important contribution. Figure IF. 141 shows the contribution 
of all parts to the final result in a logarithmic scale while Fig. IF. 151 shows 
the same in a linear scale. In all cases, i.e. low and high electron energies 
and low and high ratios between hw and Ekin,i, the main contribution comes 
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Figure F.15: Contribution of (|6T|) - (|66)) to ((60|) in a linear plot for different electron and 
photon energies {Z = 7). 
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from fl63p . It is important to state that not all contributions can be seen in 
figure IF. 141 because some of the terms have negative values which, however, 
are shown in figure IF. 151 So one might think that for Ekin,i = 100 keV and 
hw = 1 keV, equation fl66|) has the largest contribution, but as figure IF. 151 
shows, f l65|) has nearly the same absolute value, but opposite sign; therefore 
they cancel. Thus the third integral flB5]) is the most important one. The 
same holds for other electron energies and ratios between hw and Ekin,i- We 
conclude that (163|) is the dominant contribution for all relevant parameter 
values. 

Appendix G. Conservation of energy and momentum 

One can also gain information on the scattering angle Gj for high electron 
energies from the conservation of energy and momentum, 

Ei + Eg = Ef + ho, (G.l) 

p* + q = Pf + hk (G.2) 

where Eij and Pij are the energy and the momentum of the electron in the 
initial and final state, is the momentum of the photon which is related 
to its energy hw through 



h\k\ = -w, (G.3) 

c 

and Eg and q are the energy and the momentum of the virtual photon be- 
tween electron and nucleus, q changes the momentum of the nucleus. But 
the contribution to the kinetic energy can be neglected; thus Eg = and 

Ei = Ef + fko, (G.4) 
Pi-^k = P/-q. (G.5) 

Squaring flG.5p and using pj ■ k = |pi| |k| cos <(pj, k) = pi k cos 9j, the angle 
Qi is: 
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By using (lG.4p and the relativistic energy-momentum relation we get an 
expression for the momentum of the electron in the final state 



= s^jpf + fPk^ - 2nk^Jpi + mlc^ (G.7) 
which leads to 

2t\hj^ + 2qJiEi - hcuV - mld^ cos<(p/-, q) - cq^ 
cosQi = ^-^^ ^-^-^ -. G.8 

Although this is an analytical expression for the scattering angle Gj one 
should take into account that it depends on the vector q of the virtual photon 
which is not known in forehand. Thus, depending on q, only a statistical 
statement can be made about 0,. 

It is, however, possible to investigate the limit of (IG.8I) for high electron 
energies which yields 

cq 

lim cos6i = 1 + — cos<(pf, q) (G.9) 

As 6i G M <^=^ cos 6i G [—1, +1] and c, q,hco > we can conclude that 

cos<(p/,q)<0 (G.IO) 

Especially for |^ cos <(p/, q) | ^ l,6j ~ 0, i.e., the photon is mainly emitted 
in forward direction. 

If, additionally, the photon energy hw also increases more and more (for 
high electron energies) it is 

lim fl + |^cos<(p.,q)) = 1; (G.ll) 

thus 

lim Qi = 0. (G.12) 

Hence, we conclude from simple considerations about energy and momentum 
conservation that the photon is mainly scattered in forward direction if the 
energies of electron and photon are both very high. 
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Appendix H. Approximation for 

In order to obtain (11081] we calculate the derivative of f l81l) after Qf. 



d 



(Pa 



du sin QidQi 
sin^ Qi 



TT 



4-w|pj| sin6j 



sin^ 9, 



— c|pj| cosGj)^ ^ f ' ^ " ' ' Ef Ei — c\pi\cosQi 





f 2 sin Qi cos Qi{Ei — c\pi\ 


cos Qi) — 2c 


Pi 


sin^ Qi 




\ {Ei-c\pi 


cos Qi)^ 



+ 

X (4Ej + 5(e,)c2) 
+ 



2/icu; 


Pi\ 


sin^ Qi 


{Ei - c\pi 


cos 6j)^ 





2 sin 6j cos Qi{Ei 


— c\pi\ cosGj) — c\pi\ sin'^ 0j 


^/ 


{E, 


- c|pi| cos6i)2 



(H.l) 

with definition ( !77|) for 5. In order to calculate the extrema one has to set 
equation ( 1H.1I) equal to zero: 



4-a; Pi 

c 



sin 9, 



^^^^ rn n" (^i-clPilcosGi)^ 

Ef Ei — c\pi\ cos J 

5(6,) [2(EiCosei - c|p,|)(4Ej + 5(e,)c2) 

2^ca;|pi| sin^ 6i(-E'j — c|pi| cos Qi) 
2h^u^ 



+ 



(E,-c|pi|cOSe,)2^^^^ 

sin^ 6, 



+ 



Ef 



(2 cos Qi{Ei — c|pj| cos Qi 



c\pi\ sin^ Qi{Ei - c\pi\ cosGi))] . 



(H.2) 



As 6{Qi) ~ cosGj, expression flH.2p is quartic in cosG,; therefore flH.2p 
can be solved analytically in principle, but the solution will be long and 
complicated. Figure [2] also shows that the angles for maximal scattering are 
very small for relativistic electrons, therefore one can approximate cos Qi ^ 1 
and sinB, ~ G,. This leads to 



S{Q^ 



-\Pi 



-u + 2-u;|pi 
c / c 



6{Q, = 0) =: 60 



(H.3) 
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and 



7i 



= 4-a;|pi 

c 



9fi2, ,2 

e^(4Ej + 5oc^) + ^e^(E,-c|p, 



+ 60 



2{AE} + 5oc') - 2;ica;|p,|e2 + (2(Ei - c|p,|) - cIp.lG^) 



(H.4) 



with solution 



25 ohm 
E, 



{Ei - c\pi\) 



2M 



AE} + 5oc2 + ^(E, - c|p,|) - |p,|(5oc - |fc|p,|(5o 



(H.5) 



Appendix I. CH — h code 

In this section we present a C++ code which can be used to create data 
tables for the doubly differential cross sections (160|) and The code con- 
sists of three parts: In the beginning there is a header which defines all used 
physical constants and the functions for Bremsstrahlung and pair produc- 
tion. These lines have to appear in the beginning of a C++ code. Then the 
functions for Bremsstrahlung and pair production follow which can be placed 
any where else in the code. 

The function to calculate the doubly differential cross section for Bremsstrahlung 
can be called as 

brems ( E_kin , omega , Theta_i ) 

where E_kin denotes the kinetic energy of the incident electron in Joule, 
omega is the frequency of the emitted photon in s~^ and ThetaJ is the scat- 
tering angle. The function to calculate the doubly differential cross section 
for pair production can be called as 

pair_pr (E_kin , omega , Theta_p ) 

where E_kin denotes the kinetic energy of the created positron in Joule, 
omega the frequency of the incident photon in s~^ and Theta_p the angle 
between both. 

The whole code to calculate the doubly differential cross sections for 
Bremsstrahlung and pair production reads 
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// Doubly differential cross section for 
// Bremsstrahlung and pair production 

// 

// (C) Christoph Koehn, Ute Ehert 

// Centrum Wiskunde en Informatica 

//P.O. Box 94079, NL-1090 GB, Amsterdam 



#include <math . h> 

// Define physical constants 

// to be set in the beginning of the code 

#define hbar (1.054571726*pow(10, -34) ) 

// Planck's constant 
#define SpeedOfLight 299792458.0 

// Speed of light in vacuum 
#define SpeedOfLight2 pow(299792458 ,2) 

// Speed of light in vacuum, squared 
#define elcctronMass (9. 10938291*pow(10 , -31) ) 

// Mass of an electron 
#define electronCharge (1.602176565*pow(10 , -19) ) 

// Charge of an electron (without sign) 
#define rest_energy ( electronMass * SpccdOfLight2 ) 

// Rest energy of the electron 
#define alpha_finc (7.2973525698*pow(10 , -3) ) 

// fine structure constant 
#define PI 3.14159265358979 

// PI 



// Define functions for doubly differential 
// cross section for Bremstrahlung and 
// pair production , also to be set in the 
// b eginning of the code 



// Doubly differential cross section for Bremsstrahlung 
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// which depends on the kinetic energy of the incident 

electron , 

// frequency of the emitted photon and angle between 

hath 
// (in radian ) 

double brems (double E_kin, double omega , double Theta_i); 

// Doubly differential cross section for pair 
production 

// which depends on the kinetic energy of the generated 
positron , 

// frequency of the incident photon and angle between 

both 
// (in radian ) 

double pair _pr (double E_kin, double omega , double Theta_p 

); 

// 

// functions for calculating doubly differential 
// cross sections for Bremsstrahlung and pair 
// pro duction 

// Bremsstrahlung 

double brems(double E_kin , double omega , double Theta_i) 
{ 

int i ; 

// define physical quantities : energies and momenta 
double E_i=E_kin+rest -energy ; // total energy of 

incident electron 
double E_f=E_i— hbar *omega ; // total energy of 

outgoing electron 
double p_i=sqrt (pow(E_i ,2) *pow( SpeedOf Light ,— 2)— pow( 

electronMass , 2 ) * SpeedOfLight2 ) ; // momentum of 

incident electron 
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double p_f=sqrt (pow( E_f , 2 ) *pow( SpeedO flight , — 2)— pow( 
electronMass , 2 ) * SpeedOfLight2 ) ; // momentum of 
outgoing electron 

// define help functions 

int Z=7; // atomic number, can be changed 

arbitrarily 

double A=pow(Z , 2 ) *pow( alpha_fine ,3) / (4*pow( PI , 2 ) ) * p_f 

/ p _i *pow( hbar *2*PI , 2 ) /omega ; // overall 

prefactor , defined in (22) 
double B=pow(Z , 2 ) *pow( alpha_fine ,3) /(2*PI) * p_f / p_i * 

pow ( hbar *2*PI , 2 ) /omega ; // overall 

prefactor for Theta-i=0 
// delta functions defined in (51) and (52) 
double tri_l=— pow( p_i ,2)— pow( p_f ,2)— pow(hbar*pow( 

SpeedOfLight , — 1) *omega , 2 ) +2*hbar*pow( SpeedOfLight 

, — 1) *omega* p_i *cos ( Theta_i ) ; 
double t r i _2 =—2*hbar*pow( SpeedOfLight , — 1) *omega* p _f 

+2* p_i * p_f * cos ( Thet a_i ) ; 

double I[8] = {0}; // parts of the cross section 

double cross=0; // total result of cross 

section 

// all terms which are needed for the cross section 
if (Theta_i! = 0) // calculate cross section for 

all angles, except for degree 

{ 

I [0] = 2*PI*A/ sqrt (pow ( t r i_2 , 2 ) +pow(2* p_i*p_f*sin ( 
Theta_i) ,2) ) *log ( (pow( tri_2 , 2 ) +pow(2* p _i * p_f * sin 
(Theta_i) ,2) -( t ri _l+t r i _2 ) * sqrt (pow( tri_2 ,2)+pow 
(2*p_i*p_f*sin( Theta_i ) ,2)) + tri_l*tri_2) /(-pow( 
tri_2 , 2 )— pow(2* p _i * p_f * sin ( Thet a_i ) ,2)— (tri_l — 
tri_2 )*sqrt (pow( t r i_2 , 2 )+pow(2* p_i*p_f*sin ( 
Theta_i) ,2) ) + tri_l *tri_2 ) ) *(l+SpeedOfLight*tri_2 
/( p_f *( E_i — p_i * SpeedOfLight * cos ( Theta_i ) ) )— pow( 
p_i* SpeedOfLight* sin ( Thcta_i) / ( E_i-SpeedOfLight * 
p_i*cos(Theta_i)) ,2) — 2*pow( hbar *omega ,2) *p_f * 
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tri_2 /( Speed Of Light *(E_i— SpeedOf Light *p_i* cos ( 
Thet a_i ) ) * (pow ( t ri _2 , 2 )+pow(2* p_i*p_f*sin ( 
Theta_i) ,2) ) ) ) ; 

I[l] = -2*PI*A*SpeedOfLight/(p_f *(E_i-SpeedOfLight* 
p_i*cos (Theta_i) ) )*log (( E_f+SpeedOfLight*p_f ) /( 
E_f-SpeedOfLight*p_f ) ) ; 

I [2]=2*PI*A/sqrt (pow ( t r i _2 * E_f+t r i _ 1 * p _f * 

SpeedOf Light , 2 ) +pow(2* rest_energy*p_i*p_f*sin ( 
Theta_i) ,2) ) *log ( ( ( E_f+p_f *SpeedOfLight ) *(pow(2* 
p_i*p_f*sin (Theta_i) ,2) * ( E_f-p _f * SpeedOfLight ) +( 
tri_l + tri_2) *(( t ri _2 * E_f+tri _1 * p_f* SpeedOfLight ) 
— sqrt (pow( tri_2 *E_f+tri_l * p_f * SpeedOfLight ,2) + 
pow (2* rest_energy*p_i*p_f*sin (Theta_i) ,2) ) ) ) ) / (( 
E_f—p_f * SpeedOfLight ) * (pow (2* p_i*p_f*sin (Theta_i 
) ,2)*(-E_f-p_f*SpeedOfLight)+(tri_l-tri_2 ) *(( 
tri_2 *E_f+tri_l *p_f* SpeedOfLight ) — sqrt (pow( tri_2 
*E_f+tri_l * p_f * SpeedOfLight ,2)+pow(2* rest_energy 
*p_i*p_f*sin (Theta_i) ,2 ))))))*(-( (pow( t ri _2 ,2) + 
pow (2* p_i*p_f*sin(Theta_i) ,2))* (pow( E_f ,3 )+E_f * 
pow( p_f* SpeedOfLight ,2) )+p _f * SpeedOfLight *(2*( 
pow(tri_l , 2 )— pow(2* p_i * p_f * sin ( Theta_i ) ,2))*E_f* 
p_f * SpeedOf Light + t r i _ 1 * tri_2 *(3*pow( E_f ,2)+pow( 
p_f*SpeedOfLight ,2) ))) /(pow( tri_2 *E_f+tr i_l *p_f * 
SpeedOfLight , 2 ) +pow(2* rest_energy*p_i*p_f*sin ( 
Theta_i) ,2)) -SpeedOfLight *( tri_2 *E_f+tri_l *p_f* 
SpeedOfLight ) /( p_f * ( E_i— SpeedOfLight * p_i * cos ( 
Thcta_i))) -(pow(2* E_i*p_f ,2) *(2*pow( tri_2 *E_f+ 
tri_l*p_f* SpeedOfLight , 2 )— pow(2* electronMass* 
SpeedOfLight 2 * p_i*p_f*sin (Theta_i) ,2) ) *(tri_l* 
E_f+tri_2*p_f* SpeedOfLight ) ) / (pow(pow( tri_2 *E_f+ 
tri_l *p_f* SpeedOfLight ,2)+pow(2* rest_energy*p_i* 
p_f *sin ( Theta_i ) ,2) ,2) ) +(8*pow( p_i* p_f * 
rest_energy*sin(Theta_i) ,2) * (pow( E_i , 2 ) +pow( E_f 
, 2 ) ) — 2*pow ( hbar* omega* p_i*sin(Theta_i) ,2)*p_f* 
SpeedOfLight *(tri_2 *E_f+tri_l *p_f* SpeedOf Light ) 
+2*pow( hbar*omega , 2 ) * p_f *pow( electronMass ,2) *pow 
(SpeedOfLight ,3) * ( t r i _2 * E_f+t r i _ 1 * p_f * 
SpeedOfLight ) ) / ( ( E_i -SpeedOfLight * p_i *cos ( 
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Theta_i) ) *(pow( tri_2 *E_f+tri_l * p_f*SpeedO flight 
, 2 )+pow(2* rest_energy*p_i*p_f*sin (Theta_i) ,2) ) ) ) 

I[3] = -4*PI*A*p_f*SpeedOfLight *( tri_2*E_f+tri_l *p_f* 
SpeedOf Light ) /(pow( tri_2 *E_f+tri_l *p_f * 
SpeedOf Light , 2 ) +4*pow( rest_energy*p_i*p_f*sin( 
Theta_i) ,2) ) -16*PI*pow( E_i ,2) *pow(p_f ,2) *A*pow( 
t r i _2 * E_f+t r i _ 1 * p _f * SpccdOfLight ,2) /pow(pow( 
tri_2 *E_f+tri_l *p_f*SpeedOfLight ,2) +4*pow( 
rest -energy * p_i * p_f* sin ( Theta_i ) ,2) ,2) ; 

I [4]=4*PI*A/((-pow( tri_2 ,2)+pow( tri_l , 2 )-pow(2* p_i * 
p_f*sin(Theta_i) ,2)) * (pow ( t r i _2 * E_f+t r i _ 1 * p _f * 
SpeedOf Light , 2 ) +pow(2* rest_energy*p_i*p_f*sin ( 
Theta_i) ,2)))* (pow( hbar *omega* p_f , 2 ) / ( E.i- 
SpeedOfLight * p_i * cos (Theta_i) ) * ( E_f * (2*pow( tri_2 
,2) *(pow( tri_2 ,2)— pow( tri_l , 2 ) ) +8*pow( p_i * p_f * 
sin(Theta_i) ,2) *(pow( tri_2 ,2)+pow( tri_l ,2)))+p_f 
*SpeedOfLight*(2*tri_l *tri_2 *(pow( tri_2 ,2)— pow( 
tri_l ,2)) — tri_l*tri_2 *pow(4* p_i*p_f*sin (Theta_i) 
,2) ) ) / (pow( tri_2 ,2)+pow(2* p_i*p_f*sin( Theta_i ) 
,2) )+2*pow(hbar*omega*p_i*sin (Theta_i) ,2) *(2* 
tri_l*tri_2 * p _f* Speed Of Light +2*pow( tri_2 ,2)*E_f 
+8*pow( p_i * p_f * sin ( Theta_i ) ,2)*E_f) /(E_i— p_i* 
SpeedOfLight*cos (Theta_i) )+2*pow( E_i*p_f ,2) *(2*( 
pow(tri_2 ,2)— pow( tri_l ,2) ) *pow( tri_2 *E_f+tri_l * 
p_f *SpeedOfLight ,2)+8*pow( p_i*p_f*sin (Theta_i) 
,2) *((pow( tri_l ,2)+pow( tri_2 , 2 ) ) * (pow( E_f , 2 )+pow 
( p_f *SpeedOfLight ,2) ) +4* t r i _ 1 * t r i _2 * E_f * p_f * 
Speed Of Light ) ) / (pow( tri_2 *E_f+tri_l *p_f * 
SpeedOfLight ,2)+pow(2*rest_energy*p_i*p_f*sin ( 
Theta_i) , 2) ) +8*pow( p_i * p_f * sin ( Theta_i ) ,2)*(pow( 
E_i ,2)+pow(E_f ,2) ) *( tri_2 *p_f* SpeedOf Light + tri_l 
*E_f) /(E_i-SpeedOfLight*p_i*cos (Theta_i) ) ) ; 

I [5] = 1 6* PI*pow( E_f , 2 ) *pow( p_i , 2 ) *pow( sin (Theta_i) 
,2) *A/(pow(E_i-SpeedOfLight*p_i*cos(Theta_i) ,2) 
*(— pow(tri_2 ,2)+pow( tri_l , 2 ) — 4*pow( p_i ,2)*pow( 
p_f ,2)*pow(sin(Theta_i) ,2))) ; 
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else // calculate cross section for degree 

{ 

I[0] = -2*B*p_f*SpeedOfLight /( tri_2 *E_f+tri_l *p_f * 
SpeedOfLight) ; 

I [1] =B*pow( p_f *SpeedOfLight ,2) *(— pow( tri_l ,2)+pow( 
tri_2 ,2) ) / ( tri_2 *pow( t r i _2 * E_f+t r i _ 1 * p_f * 
SpeedOfLight , 2) ) *log ( ( t ri _l+t r i _2 ) / ( tri_l -t r i _2 ) 

); 

I [2]=B*(2* tri_l *E_f*p_f *SpeedOfLight+tri_2 *(pow(E_f 
,2)+pow(p_f*SpeedOfLight ,2) ) ) /(pow( tri_2 *E_f+ 
tri_l*p_f*SpeedOfLight ,2) ) * log ( ( E_f+p_f * 
SpeedOfLight)/(E_f-p_f*SpeedOfLight)) ; 

I [3] = -B*pow((4* E_i*p_f ) / ( tri_2 *E_f+tri_l *p_f * 
SpeedOfLight) ,2) ; 

I [4] = -8*B*pow( E_i*p_f ,2) *( t r i _ 1 * E_f+t r i _2 * p_f * 
SpeedOfLight ) /(pow( tri_2 *E_f+tri_l *p_f * 
SpeedOfLight , 3 ) ) *log ( ( ( tri_l -t ri _2 ) *(E_f-p_f* 
SpeedOfLight) ) / (( tri_l + tri_2 ) *(E_f+p_f* 
SpeedOfLight) ) ) ; 

1 [5] = — B*pow(2* hbar * p_f *omega , 2 ) / ( ( t r i _2 * E_f+t r i _1 * 
p_f *SpeedOfLight ) *(E_i-p_i*SpeedOfLight ) *tri_2 ) ; 

I [6] = 2*B*pow(hbar*p_f *omega ,2)*(2*tri_l*tri_2 *E_f+ 
p_f*SpeedOfLight *(pow( tri_l , 2 )+pow( tr i_2 ,2)))/( 
pow( t r i _2 *E_f+t r i _1 * p _f * SpeedOfLight , 2 ) * ( E_i — p_i 
*SpeedOfLight)*pow(tri_2 ,2) ) *log ( ( tri_l+tri_2 ) / ( 
tri_l -tri_2 ) ) ; 

I [7] = 2*B*pow ( hbar * omega , 2 ) * p _f * (pow ( E_f . 2 )— pow( p_f * 
SpeedOfLight ,2) ) / (pow( t r i _2 * E_f+t r i _ 1 * p _f * 
SpeedOfLight ,2) *(E_i-SpeedOfLight*p_i )* 
SpeedOfLight ) *log ( ( E_f-p_f *SpeedOfLight ) /( E_f+ 
p_f*SpeedOfLight)) ; 

} 

// sum up all terms to get value for cross section 
for (i=0;i<=7;i++) 

{ 

cross=cross+I [ i ] ; 

} 
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return cross ; // give back result of cross section 
calculation 

} 

// pair pro due ti on 

double pair_pr (double E_kin, double omega , double Theta_p 
) 

{ 

int i ; 

// define physical quantities: energies and momenta 
double E_p=E_kin+rest -energy ; // total energy 

of created positron 
double E_m=(hbar*omega— rest -energy )—E_kin ; // 

total energy of created electron 
double p_p=sqrt (pow( E_p , 2 ) *pow( SpeedOf Light , — 2)— pow( 

electronMass , 2 ) * SpeedOfLight2 ) ; // momentum of 

created positron 
double p_m=sqrt (pow(E_m,2) *pow( SpeedOfLight ,— 2)— pow( 

electronMass , 2 )* SpeedOfLight2 ) ; // momentum of 

created electron 

// define help functions 

int Z=7; // atomic number, can be changed 

arbitrarily 

double A=pow(Z , 2 ) *pow( alpha_f ine , 3 ) * SpeedOfLight2 *p_p 

*p_m/ (4*pow(PI ,2) *(hbar*2*PI) *pow( omega ,3) ) ; 

// overall prefactor , defined in (107) 
double B=pow(Z , 2 ) *pow( alpha_fine ,3) *SpeedOfLight2*p_p 

*p_m / (2*PI*(hbar*2*PI) *pow( omega , 3 ) ) ; 

// overall prefactor for Theta^p—O and Theta-p—180 
degree 

// delta functions defined in (108) and (109) 
double tri_l=— pow(p_p , 2 )— pow(p_m , 2 )— pow( hbar*pow( 

SpeedOfLight , — 1) *omega , 2 ) +2*hbar*pow( SpeedOfLight 

, — 1) *omega*p_p*cos ( Theta_p ) ; 
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double tri_2=2*hbar*pow(SpeedOfLight , — 1) *omega*p_m— 2* 
p_p*p_m*cos (Theta_p) ; 

double I[8] = {0}; // parts of the cross section 

double cross=0; // total result of cross 

section 

// all terms which are needed for the cross section 
if (Theta_p! = &fe Theta_p! = PI) // calculate 

cross section for all angles, except for and 180 
degree 

{ 

I [0]=2*PI*A/ sqrt (pow ( t r i _2 , 2 ) +pow(2* p_p*p_m* sin ( 
Tlieta_p) ,2) ) * log ( (pow( tri_2 ,2)+pow(2*p_p*p_m*sin 
(Theta_p) ,2 ) -( t r i _l+t r i _2 ) * sqrt (pow ( t ri _2 ,2)+pow 
(2* p_p*p_m* sin (Theta_p) ,2))+tri_l*tri_2) /(— pow( 
tri_2 , 2 )— pow(2* p_p*p_m* sin ( Theta_p ) ,2)— (tri_l — 
tri_2 )*sqrt (pow( t r i_2 , 2 )+pow(2* p_p*p_m* sin ( 
Theta_p) ,2) ) + tri_l * tri_2 ) ) *(-l -SpeedOfLight * 
tri_2 / (p_m* (E_p—p_p* Speed Of Light * cos (Theta_p) ) ) + 
pow(p_p*SpeedOfLight*sin (Theta_p) / (E_p— 
SpeedOfLight *p_p* cos (Theta_p) ) ,2) — 2*pow(hbar* 
omega , 2 ) *p_m* tri_2 / (SpeedOfLight* (E_p— 
SpeedOfLight *p_p* cos ( Theta_p ) ) * (pow( t r i_2 , 2 ) +pow 
(2* p_p*p_m* sin (Theta_p) ,2) ) ) ) ; 

I [l]=2*PI*A*SpeedOfLight /(p_m*(E_p-SpeedOfLight*p_p 
*cos (Theta_p) ) )*log ( (E_m+SpeedOf Light *p_m) /(E_m- 
SpeedOfLight *p_m) ) ; 

I[2]=2*PI*A/sqrt (pow ( t r i _2 *E_ni+t r i _ 1 *p_m* 

SpeedOfLight , 2 ) +pow(2* rest_energy *p_p*p_m* sin ( 
Theta_p ) ,2) ) * log ( (( E_m+p-m* SpeedOfLight ) *(pow(2* 
p_p*p_m* sin(Theta_p) ,2)* (E_ni—p_m* SpeedOfLight ) +( 
tri_l + tri_2 ) *(( tri_2 *E_ni+t r i _1 *p_m* SpeedOfLight ) 
— sqrt (pow( t r i _2 *E_ni+t r i _ 1 *p_m* SpeedOfLight , 2 ) + 
pow (2* rest -energy *p_p*p_m* sin (Theta_p) ,2) ) ) ) ) / (( 
E_m.—p_m* SpeedOfLight ) * (pow(2* p_p*p_m* sin ( Theta_p 
) ,2) *(-E_m-p_m* SpeedOfLight )+(tri_l - tri_2 ) *(( 
t r i _2 *E_ni+t r i _ 1 *p_ni* SpeedOfLight ) — sqrt (pow( t r i _2 
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*E_m+tri_l *p_m*SpeedOfLiglit , 2)+pow(2* rest_energy 
*p_p*p_m* sin ( Theta_p ) ,2) )))))*(( (pow( tri_2 ,2) + 
pow(2* p_p*p_m* sin ( Theta_p ) , 2) ) * (pow (E_m, 3 ) +E_iii* 
pow(p_m*SpeedOfLight ,2) )+p_m*SpeedOfLight *(2*( 
pow(tri_l ,2)— pow(2*p_p*p_m*sin (Theta_p) ,2))*E_m* 
p_m*SpeedOfLight + tri_l * tri_2 *(3*pow(E_m,2)+pow( 
p_m* SpeedOf Light , 2 ) ) ) ) / (pow( t r i _2 *E_m+t r i _1 *p_m* 
SpecdOf Light , 2 ) +pow(2>K rest_energy *p_p*p_m* sin ( 
Theta_p ) ,2) )+SpeedOf Light *( tri_2 *E_m+tri_l *p_m* 
SpeedOfLight ) / (p_m*(E_p— SpeedOfLight*p_p*cos ( 
Theta_p) ) )+(pow(2*E_p*p_m,2) *(2*pow( tri_2 *E_m+ 
tri_l *p_m* SpeedOfLight ,2)— pow (2* electronMass* 
SpeedOfLight 2 *p_p*p_m* sin(Theta_p) ,2))*(tri_l* 
E_m+tri_2 *p_m*SpeedOfLight ) ) / (pow(pow( tri_2 *E_m.+ 
tri_l *p_m* SpeedOfLight ,2)+pow(2* rest_energy*p_p* 
p_m*sin (Theta_p) ,2) ,2) )+(— 8*pow(p_p*p_m* 
rest_energy*sin(Thcta_p) ,2) * (pow( E_p , 2 ) +pow(E_m 
, 2 ) ) — 2*pow ( hb ar* omega* p_p* sin ( Theta_p ) ,2 ) *p_m* 
SpeedOfLight *( tri_2 *E_m+tri_l *p_m*SpeedOfLight ) 
+2*pow( hbar*omega , 2 ) *p_m*pow( electronMass , 2 ) *pow 
(SpeedOfLight ,3) *( tri_2 *E_m+tri_l *p_m* 
SpeedOfLight ) ) / (( E_p— SpeedOfLight *p_p * cos ( 
Theta_p ) ) * (pow ( t r i _2 *E_ni+t r i _1 *p_m* SpeedOfLight 
, 2 ) +pow(2* rest - energy *p_p*p_m* sin ( Theta_p ) ,2 ) ) ) ) 

I [3] =4* PI *A*p_m* SpeedOfLight *( tri_2 *E_m+tri_l *p_m* 
SpeedOfLight ) / (pow( tr i_2 *E_m+tri_l *p_m* 
SpeedOfLight , 2 ) +4*pow( rest_energy *p_p*p_m* sin ( 
Theta_p ) , 2) ) +16*PI*pow( E_p , 2 ) *pow(p_m , 2 ) *A*pow( 
tri_2 *E_m-|-tri_l *p_m*SpeedOfLight ,2) /pow(pow( 
t r i _ 2 *E_m+t r i _ 1 *p_m* SpeedOfLight , 2 ) +4*pow ( 
rest _energy *p_p*p_m* sin ( Theta_p ) ,2) ,2) ; 

I [4] =4*PI*A/((-pow( tri_2 ,2)+pow( tri_l ,2)-pow(2*p_p* 
p_m* sin ( Theta_p ) ,2 ) ) * (pow ( t r i _2 *E_m+t r i _ 1 *p_m* 
SpeedOfLight , 2 ) +pow(2* rest_energy *p_p*p_m* sin ( 
Theta_p) ,2) ) ) * (pow( hbar*omega*p_m , 2 ) / (E_p— 
SpeedOfLight *p_p*cos (Thcta_p) ) * (E_m* (2*pow( tri_2 
,2) *(pow( tri_2 , 2 )— pow( tr i _1 ,2) )+8*pow(p_p*p_m* 
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sin(Theta_p) ,2) *(pow(tri_2 ,2)+pow( tri_l ,2) ) )+P-m 
*SpeedOfLight*(2* tri_l *tri_2 * (pow( t r i _2 , 2 )— pow( 
tri_l ,2)) — tri_l*tri_2 *pow(4* p_p*p_m>K sin ( Theta_p ) 
,2 ) ) ) / (pow ( t ri _2 , 2 )+pow(2* p_p *p_m* sin ( Theta_p ) 
,2) )+2*pow(hbar*omega*p_p* sin (Theta_p) ,2) *(2* 
tri_l*tri_2 *p_m* SpeedOf Light +2*pow( t ri _2 , 2 ) *E_m 
+8*pow( p_p*p_m* sin(Theta_p) ,2) *E_m) / ( E_p— p_p* 
Speed Of Light *cos(Theta_p)) — 2*pow( E_p*p_m , 2 ) * (2* ( 
pow(tri_2 , 2 )— pow( tri_l ,2) ) *pow( tri_2 *E_m+tri_l * 
p_m* SpeedOf Light ,2) +8*pow( p_p*p_m* sin (Theta_p) 
, 2 ) * ( ( pow ( t r i _ 1 , 2 ) +pow (tri_2 ,2))*( pow ( E_ni , 2 ) +pow 
(p_m* SpeedOf Light , 2 ) ) +4* tri_l*tri_2 *E_ni*p_ni* 
Speed Of L ight ) ) / (pow( tri_2 *E_m+tri_l *p_m* 
SpeedOf Light , 2 ) +pow(2* rest_energy *p_p*p_m* sin ( 
Theta_p) ,2) )— 8*pow(p_p*p_m*sin (Theta_p) ,2)*(pow( 
E_p ,2)+pow(E_m,2) ) *( tri_2 *p_m*SpeedOfLight+tri_l 
*E_ni) / (E_p— SpeedOfLight*p_p*cos (Theta_p) ) ) ; 
I [5] = — 16* PI*pow(E_m, 2 ) *pow( p_p , 2 ) *pow ( sin ( Theta_p ) 
,2 ) >kA/ (pow (E_p— Speed Of Light *p_p*cos(Theta_p) ,2) 
*(— pow(tri_2 ,2)+pow( tri_l ,2) — 4*pow(p_p ,2) *pow( 
p_m,2) *pow( sin (Theta_p) ,2) ) ) ; 

} 

else 
{ _ 

if ( Theta_p! = PI) // calculate cross section for 
degree 

{ 

I [0] =2>KB*p_m* SpeedOfLight / ( tri_2 *E_ni+tri_l *p_m* 

SpeedOfLight) ; 
I [1] = — B*pow(p_m*SpeedOfLight ,2) *(— pow( tri_l ,2) + 

pow( tri_2 ,2))/(tri_2 *pow( tri_2 *E_ni+tri_l *p_m* 

SpeedOfLight , 2 ) ) * log ( ( t r i _ 1+t r i _2 ) / ( t ri _1 - 

tri-2)); 

I [2] = — B*(2* tri_l *E_m*p_m* SpeedOf Light+t r i _2 *(pow( 
E_ni, 2 )+pow(p_m* SpeedOfLight ,2) ) ) / (pow( tri_2 * 
E_ni+tri_l *p_ni* SpeedOfLight ,2) ) *log ( (E_ni+p-m* 
SpeedOfLight ) /(E_m-p_m* SpeedOfLight ) ) ; 
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I [3] =B*pow ( (4* E_p*p_m) / ( t r i _2 *E_m+t r i _ 1 *p_m* 

SpeedOfLight) ,2) ; 
I [4] = 8*B*pow( E_p*p_m , 2 ) * ( t r i _ 1 *E_m+t r i _2 *p_m* 

SpeedOfLight ) / (pow( tri_2 *E_m+tri_l *p_m* 

SpeedOfLight ,3) ) *log ( ( ( tri_l -tri_2 ) *(E_m-p_m* 

SpeedOfLight ) ) / ( ( t r i _ 1+t r i _2 ) *(E_m+p_m* 

SpeedOfLight))) ; 
I [5] = — B*pow(2*hbar*p_m*omcga ,2) / ( ( tri_2 *E_m+tri_l 

*p_m* SpeedOfLight ) *(E_p— p_p*SpeedOfLight ) * 

tri-2) ; 

I [6] = 2*B*pow( hbar*p_m*omega ,2)*(2*tri_l*tri_2 *E_m 
+P-m*SpeedOfLight *(pow( tri_l ,2)+pow( tri_2 ,2))) 
/ (pow( t r i _2 *E_m+t r i _ 1 *p_m* SpeedOfLight , 2 ) * (E_p 
—p_p* SpeedOfLight ) *pow( tri_2 ,2))*log((tri_l + 
tri_2 ) /( tri_l -tri_2 ) ) ; 

I [7] = 2*B*pow(hbar*omega ,2) *p_m*(pow(E_m,2)— pow( 
p_m* SpeedOfLight ,2) ) / (pow( tri_2 *E_m+tri_l *p_m* 
SpeedOfLight ,2) * (E_p-SpeedOfLight *p_p ) * 
SpeedOfLight ) * log ( (Ej2i—p_m* SpeedOfLight ) / (E_m.+ 
p_m* SpeedOfLight ) ) ; 

} 

else // calculate cross section for 180 

degree 

{ 

I [0] = 2 *B*p_m* SpeedOfLight / ( t r i _2 *E_m-|-t r i _ 1 *p_m* 

SpeedOfLight) ; 
I [l] = —B*pow(p_m* SpeedOfLight ,2) *(— pow( tri_l ,2) + 

pow( tri_2 ,2))/(tri_2 >Kpow( tr i_2 *E_m+tr i_l *p_m* 

SpeedOfLight ,2) ) *log ( ( tri_l + tri_2 ) /( tri_l - 

tri-2)); 

I [2] = — B*(2* tri_l *E_m*p_m*SpeedOfLight+tri_2 *(pow( 

E_m , 2 ) +pow (p_m* SpeedOfLight , 2 ) ) ) / (pow( t r i _2 * 
Ejii+t r i _ 1 *p_m* SpeedOfLight , 2 ) ) * log ( (Ejn+P-m* 
SpeedOfLight) / (E_m— p_m* SpeedOfLight)) ; 

I [3] =B*pow ( (4* E_p*p_m) / ( t r i _2 *E_m-|-t r i _ 1 *p_m* 
SpeedOfLight) ,2) ; 

I [4] = 8*B*pow( E_p*p_m , 2 ) * ( t r i _ 1 *E_m+t r i _2 *p_in* 
SpeedOfLight ) / (pow( tri_2 *E_m+tri_l *p_m* 
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SpeedOfLight ,3) ) * log ( ( ( tri_l -tri_2 ) *(E_m-p_m* 
SpccdOfLight ) ) /(( tri_l+tri_2 ) *(E_m+p-m* 
SpeedOfLight))) ; 
I [5] = — B*pow(2*hbar*p_m*omega ,2) / ( ( tri_2 *E_m.+ tri_l 
*p_m* SpeedOfLight ) *( E_p+p_p* SpeedOfLight ) * 
tri_2) ; 

I [6] = 2*B*pow ( hbar*p_m* omega ,2)*(2*tri_l*tri_2 *E_m 
+p_m*SpeedOfLight *(pow( tri_l ,2)+pow( tri_2 ,2))) 
/ (pow( t r i _2 *E_m+t r i _ 1 *p_m* SpeedOfLight , 2 ) * ( E_p 
+p_p* SpeedOfLight ) *pow( tri_2 ,2))*log((tri_l + 
tri_2 ) /(tri_l -tri_2 ) ) ; 

I [7] = 2*B*pow ( hbar* omega , 2 ) *p_m* (pow(E_m, 2 )— pow( 
p_m* SpeedOfLight ,2) ) / (pow( tri_2 *E_m+tri_l *p_m* 
SpeedOfLight ,2) * ( E_p+SpeedOfLight *p_p ) * 
SpeedOfLight ) * log ( (E_m—p_m* SpeedOfLight ) / (E_m+ 
p_m*SpeedOfLight ) ) ; 

} 

} 



// sum up all terms to get value for cross section 
for (i=0;i<=7;i++) 

{ 

cross = cross+I [ i ] ; 

} 

return cross ; // give hack result of cross section 
calculation 

} 
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